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The  objective  of  the  present  program  is  to  study  the  structure  and 
response  of  steady  and  unsteady  laminar  premixed  and  nonpremixed 
flames  in  reduced  and  elevated  pressure  environments  through  (a) 
ron-intrusi ve  experimentation,  (b)  computational  simulation  using 
detailed  flame  and  kinetic  codes,  and  (c)  asymptotic  analysis  with 
reduced  kinetic  mechanisms.  During  the  reporting  period  progress 
has  been  made  in  the  following  projects:  (1)  a  theoretical  and 
experimental  study  of  unsteady  diffusion  flames;  (2)  a 
computational  and  experimental  study  of  hydroge  /air  diffusion 
flames  at-  sub-  and  super-atmospheric  pressures;  (3)  an  asymptotic 
analysis  of  the  structure  of  premixed  flames  with  volumetric  heat 
loss;  (4)  asymptotic  analyses  of  ignition  in  the  supersonic 
hydrogen/air  mixing  layer  with  reduced  mechanisms;  (5)  a  new 
numerical  algorithm  for  generating  the  ignition-extinction  S- 
curves .  A  total  of  three  reprints  are  appended. 
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1 .  Introduction 

The  objective  of  the  present  program  is  to  study  the  structure  and 
response  of  steady  and  unsteady  laminar  premixed  and  nonpremixed 
flames  in  reduced  and  elevated  pressure  environments  through  (a) 
non-intrusive  experimentation,  (b)  computational  simulation  using 
detailed  flame  and  kinetic  codes,  and  (c)  asymptotic  analysis  with 
reduced  kinetic  mechanisms.  During  the  reporting  period  progress 
has  been  made  in  the  following  projects;  (1)  a  theoretifcal  and 
experimental  study  of  unsteady  diffusion  flames;  (2)  a 
computational  and  experimental  study  of  hydrogen/air  diffusion 
flames  at  sub-  and  super-atmospheric  pressures;  (3)  an  asymptotic 
analysis  of  the  structure  of  premixed  flames  with  volumetric  heat 
loss;  (4)  asymptotic  analyses  of  ignition  in  the  supersonic 
hydrogen/air  mixing  layer  with  reduced  mechanisms;  (5)  a  new 
numerical  algorithm  for  generating  the  ignition-extinction  S- 
curves .  These  accomplishements  are  briefly  discused  in  the 
following . 


2 .  Studies  on  Unsteady  Diffusion  Flames 

A  crucial  influence  on  the  flame  behavior  which  so  far  has  not 
been  adequately  addressed  is  the  effect  of  unsteadiness  of  the 
environment  on  the  flame  behavior.  This  issue  is  of  particular 
relevance  to  the  modeling  of  turbulent  flames  through  the  concept 
of  laminar  flamelets.  These  flamelets  are  subjected  to 
fluctuating  flows  with  various  intensities  of  straining,  and  it  is 
reasonable  to  expect  that  the  flame  would  respond  differently  in 
an  oscillating  strained  flow  field  than  in  a  steady  strained  flow 
field. 

During  the  reporting  period  we  have  analyzed  the  response  of 
a  counterflow  diffusion  flame  subjected  to  an  oscillating  strain 
rate,  using  large  activation  energy  asymptotics.  The 

characteristic  oscillation  time  of  practical  interest  is  found  to 
be  of  the  same  order  as  the  characteristic  diffusion  time  of  the 
flame,  so  that  the  flame  structure  consists  of  a  quasi-steady 
reactive-diffusive  layer  embedded  in  the  outer  unsteady- 
convective-diffusive  zones.  A  linear  analysis  is  conducted  by 
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assuming  that  the  amplitude  of  the  strain  rate  oscillation  is 
small  relative  to  the  mean  strain  rate. 


Figure  1  shows  the  real  parts  of  the  fluctuations  in  heat 
release  as  a  function  of  the  frequency  of  oscillation,  when  the 
mean  flame  is  either  near  equilibrium  or  near  extinction.  It  has 
been  found  that  the  flame  response  is  controlled  mainly  by  two 
effects:  (a)  tne  response  of  the  convective  mass  flux  into  the 
reaction  sheet,  which  is  directly  related  to  the  flow-field 
variation  applied  at  the  boundary,  and  (b)  the  response  of  the 
reaction  sheet  to  adjust  the  reduced  residence  time  due  to  finite- 
rate  chemistry.  For  flames  near  equilibrium,  the  former  effect 
tends  to  be  dominant,  so  that  the  response  of  the  net  heat  release 
is  in  phase  (i.e.  positive  real)  with  the  strain  rate  oscillation. 
For  flames  near  extinction,  however,  the  finite-rate  chemistry 
effect  overtakes  the  fluid-dynamic  effect  such  that  increasing 
strain  rate  leads  to  a  reduction  of  the  reactivity  of  the  flame 
during  the  oscillatory  cycle.  As  such,  the  net  heat  release 
response  of  the  near-extinction  flame  becomes  out  of  phase  with 
the  strain  rate  oscillation  in  the  sense  of  the  Rayleigh's 
criterion.  Results  of  the  present  study  suggest  the  possibility 
that  the  unsteady  characteristics  of  the  near-extinction  diffusion 
flame  can  be  significantly  different  from  those  in  the  Burke- 
Schumann  limit. 

A  parallel  experimental  study  is  also  in  progress.  Figure  2 
shows  the  schematic  of  the  experimental  setup.  To  achieve  the 
instantaneous  measurement  of  the  unsteady  counterflow  flame,  the 
laser  beam  is  modulated  by  an  optical  chopper  at  the  applied 
perturbation  frequency.  Then  the  signal  is  delayed  and  sent  to 
the  loudspeakers,  which  generate  the  acoustic  perturbation  of  the 
flow  in  the  nozzle.  By  comparing  the  signal  sent  to  the 
loudspeakers  with  the  chopped  laser  signal  detected  by  the 
photodiode,  the  entire  applied  sine  wave  can  be  mapped  out.  We 
shall  measure  the  instantaneous  snapshots  of  the  strain  rate  field 
and  flame  responses  in  order  to  provide  more  deterministic 
information  of  the  unsteady  flame  phenomena. 
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The  theoretical  part  of  this  work  is  reported  in  Publication 
No.  1. 

3.  Pressure  Effects  on  Extinction  of  Hydrogen-Air 

Diffusion  Flames 

Recent  interests  in  high-speed  aero-propulsion  have  led  to 
considerable  research  on  hydrogen/oxygen  chemistry,  as  well  as  its 
coupling  to  fluid  flows.  Because  of  the  high-speed  nature  of  the 
flow,  the  available  residence  time  for  mixing  and  chemical 
reaction  is  significantly  reduced.  Furthermore,  the  combustion 
chamber  within  aero-engines  not  only  operates  at  higher  pressures, 
but  the  chamber  pressure  can  also  undergo  strong  fluctuations  such 
that  chemical  kinetics  and  the  flame  behavior  can  be  significantly 
modified  from  those  at  the  atmospheric  condition.  Consequently, 
it  is  of  importance  to  understand  the  ignition  and  extinction 
phenomena  involving  hydrogen/oxygen  mixtures  under  variable 
pressures . 

We  have  conducted  LDV  measurements  of  local  extinction  strain 
rates  of  nonpremixed  counterflow  flames  of  diluted  hydrogen 
against  air,  at  pressures  of  0.5  to  1.0  atmosphere.  The  measured 
data  compare  well  with  results  obtained  from  computational 
simulation  with  detailed  chemistry  and  transport.  We  have 
subsequently  performed  additional  computational  studies  of  the 
pressure  effect  on  flame  extinction.  Figure  3  shows  extinction 
flame  temperatures  calculated  for  a  13%  hydrogen  mixture  as  a 
function  of  the  system  pressure.  The  relationship  of  extinction 
temperatures  with  pressure  for  this  system  exhibits  the  familiar 
non-monotonic  "Z"  shaped  dependency  as  observed  for  the 
homogeneous  hydrogen/oxygen  explosion  limits.  This  implies  that 
an  increase  in  pressure  could  render  a  mixture  to  change  from  an 
extinguished  state,  to  a  burning  state,  and  back  to  an 
extinguished  state.  This  behavior  is  explained  on  the  basis  of 
the  intrinsic  chain  branching-termination  kinetics  of  hydrogen 
oxidation. 

This  work  is  reported  in  Publication  No.  2. 
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4.  Asymptotic  Analysis  of  Prsmixad  Flamas  with  Volumatric 

Hast  Loss 

A  classical  model  problem  for  the  study  of  premixed  flame 
extinction  is  that  of  Spalding,  who  analyzed  a  one-dimensional 
freely-propagating  flame  with  a  temperature-sensitive  one-step 
overall  reaction  and  radiative  heat  loss.  More  sophisticated 
analyses  have  since  been  performed  using  one-step  chemistry  and 
activation  energy  asymptotics.  These  studies  predict  that 
extinction  occurs  when  the  ratio  of  the  burning  rate  to  its 
adiabatic  value  is  reduced  to  e~i/2  0.61,  a  result  that  is 
insensitive  to  the  nature  of  the  heat  loss.  Subsequent  numerical 
computations  that  consider  detailed  transport  and  chemistry  also 
predict  extinction  when  the  burning  rate  is  reduced  to  roughly  60% 
of  its  adiabatic  value,  suggesting  that  this  value  may  represent  a 
universal  constant,  independent  of  reaction  and  loss  mechanisms. 

In  order  to  gain  a  better  understanding  of  the  role  of 
dominant  kinetic  parameters,  we  have  revisited  this  classical 
problem  and  have  performed  asymptotic  analyses  with  multi-step 
reaction  mechanisms.  We  first  employed  the  two-step  Zel'dovich- 
Lifi^n  mechanism  which  consists  of  a  branching  reaction  and  a 
competing  recombination  reaction,  thus  capturing  the  chain  nature 
of  real  flames.  The  analysis  again  predicts  the  normalized  burning 
rate  at  extinction  to  be  This  work  is  reported  in 
Publication  No.  3. 

We  have  also  made  considerable  progress  in  understanding  the 
extinction  characteristics  of  nonadiabatic  methane/air  flames  by 
considering  a  reduced  reaction  mechanism  that  has  been 
systematically  derived  from  larger  detailed  mechanisms.  Thus  the 
reduced  mechanism  contains  many  important  chemical  parameters  that 
are  dominant  in  the  chemistry  of  "real"  flames.  Extinction 
conditions  were  found  explicitly  in  terms  of  these  parameters. 
Results  show  that  the  critical  value  of  the  burning  rate  at 
extinction,  though  not  a  universal  constant,  varies  only  slightly, 
in  the  range  between  0.61  and  0.64.  In  figure  4  we  have  plotted 
the  nondimensional  burning  rate  as  a  function  of  the  loss 
parameter.  At  the  (extinction)  turning  point,  the  burning  rate  is 
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nearly  the  same  for  all  curves,  but  the  critical  value  of  the  heat 
loss  parameter  is  seen  to  be  less  than  that  for  one-step  models. 
The  corresponding  reduction  in  flame  temperature  needed  for 
extinction  is  therefore  substantially  less  than  that  predicted  by 
one-step  models . 

This  work  is  reported  in  Publication  Nos.  4  and  5. 

5 .  Ignition  in  the  Supersonic  Bydrogen/Air  Mixing  Iisyer 

The  development  of  the  scramjet  engine  for  supersonic  propulsion 
and  the  scramaccelerator  for  hypervelocity  projectile  launching 
has  renewed  interest  in  supersonic  combustion.  Fundamental 
studies  such  as  ignition  within  chemically-react ive  supersonic 
boundary  layer  flows  can  provide  better  understanding  of  the  key 
factors  in  developing  supersonic  combustors.  In  our  previous 
study  (Publication  No.  6)  with  one-step  reaction  chemistry  model, 
several  distinct  ignition  situations  were  identified,  and  it  was 
shown  that  ignition  can  be  greatly  facilitated  by  use  of  the 
kinetic  energy  of  the  high-speed  flow  through  viscous  heating.  In 
the  present  study  we  take  one  step  further  and  consider  the  more 
realistic  hydrogen/air  mixing  layer  using  reduced  mechanisms.  In 
particular,  two  distinct  reduced  mechanisms  are  identified 
depending  on  the  maximum  characteristic  ignition  temperature  (T*) 
relative  to  the  crossover  temperature  (Tc)  ,  at  which  the  rates  of 
the  crucial  H-02  chain  branching  and  termination  reactions  are 
equal.  Each  regime  is  further  subdivided  into  two  distinct  cases, 
namely  the  hot  stream  (P  >  fi)  case  and  the  viscous  heating  (P  <  4) 

case,  depending  on  the  relative  dominance  of  the  external  and 
Internal  ignition  sources,  where  P  and  are  parameters 

respectively  representing  the  extent  of  external  heating  to 
internal  viscous  heating.  These  four  cases  are  analyzed 
separately  using  asymptotic  technique,  and  the  analysis  properly 
captures  the  ignition  process  in  a  well-defined  manner. 

Figure  5  shows  the  minimum  ignition  distance  as  a  function  of 
the  flow  Mach  number.  Although  the  gross  trend  is  similar  to  that 
of  the  one-step  analysis,  it  is  clearly  shown  that  the  system 
response  is  significantly  enriched  when  realistic  chemistry  is 
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properly  taken  into  account.  For  example,  ignition  in  the  low 
temperature  regime  (I,  II)  is  controlled  by  a  large  activation 
energy  process,  so  that  the  ignition  distance  is  more  sensitive  to 
its 'characteristic  temperature  than  that  in  the  high  temperature 
regime  (III,  IV).  In  Fig.  6,  it  is  shown  that  the  ignition 
distance  varies  non-monotonically  with  the  system  pressure  in  the 
manner  of  the  well-known  hydrogen/oxygen  explosion  limits,  thereby 
further  substantiating  the  importance  of  chemical  chain  mechanisms 
in  this  class  of  chemically-reacting  boundary  layer  flows. 

This  work  is  reported  in  Publication  No.  7. 

6.  Numerical  Algorithm  for  Generating  S-Curves 

Most  combustion  phenomena  are  intrinsically  nonlinear  because  of 
the  associated  chain  mechanisms  and  Arrhenius  kinetics. 
Consequently,  theoretical  studies  of  their  steady-state  behavior 
frequently  yield  multiple  solutions  connected  by  turning  points 
when  a  system  response  is  plotted  versus  an  imposed  system 
parameter.  A  prominent  example  is  the  S-curve  representation  of 
quasi-one-dimensional  strained  flames  which  results  when  a  flame 
response  such  as  its  maximum  temperature  is  plotted  versus  the 
system  strain  rate.  Specifically,  the  lower  branch  of  the  S-curve 
represents  weakly-reactive  states,  the  middle  branch  unstable 
solutions,  and  the  upper  branch  vigorously-burning  states.  As 
such,  the  lower  and  upper  turning  points  can  be  respectively 
defined  as  the  ignition  and  extinction  states  of  the  system. 
Identification  and  description  of  these  critical  states  of 
transition  therefore  commands  both  fundamental  and  practical 
interest . 

When  numerical  solution  is  sought  for  the  S-curve  response, 
especially  with  detailed  chemistry,  difficulty  is  encounters  as 
one  attempts  to  continue  the  calculation  from,  say,  the  upper 
branch  to  the  middle  branch  by  negotiating  the  extinction  turning 
point  because  the  Jacobian  matrix  used  in  the  Newton  method 
becomes  singular  at  the  turning  point.  Consequently,  the 
arclength  continuation  method  has  been  applied  to  trace  through 
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the  ignition  and  extinction  turning  points  in  a  number  of  studies 
on  strained  flames. 

While  the  arclength  continuation  method  is  clearly  a  valuable 
approach  in  generating  the  S-curve,  special  skill  apparently  is 
still  needed  in  its  implementation  and  as  a  result  its  use  has  not 
been  wide  spread.  In  this  study  we  have  developed  a  new,  flame¬ 
controlling  continuation  method.  Using  the  counterflow  premixed 
and  diffusion  flames  as  examples,  the  method  capitalizes  on  the 
distinct  nature  of  the  profile  and  location  of  the  scalars  of  the 
flame  properties,  such  as  the  temperature  and  species 
concentrations,  in  response  to  changes  in  the  flow  strain  rate. 
Thus  instead  of  using  the  strain  rate  as  an  imposed  parameter  and 
the  scalars  as  the  flame  responses,  the  values  of  a  flame  scalar 
at  a  give  location  is  used  as  an  internal  boundary  condition  while 
the  strain  rate  becomes  the  flame  response.  The  method  appears  to 
be  fairly  simple  in  implementation  and  efficient  and  robust  in 
execution,  especially  in  negotiating  the  turning  points. 

This  work  is  reported  in  Phblication  No.  8. 
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Laminar  flame  propagation  with  volumetric  heat 
loss  and  chain  branching-termination  reactions 

B.  H.  CHAO 

Department  of  Mechanical  Engineering.  University  of  Hawaii.  Honolulu.  HI  96822.  U.S.A. 

and 

C.  K.  LAW 

Department  of  Mechanical  and  Aerospace  Engineering.  Princeton  University. 
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Abstract — The  steady  propagation  of  the  planar  laminar  premixed  flame  in  the  doubly-infinite  domain, 
with  chain  branching-termination  reactions  and  weak  volumetric  heat  loss,  is  studied  using  activation 
energy  asymptotics.  Two  flame  propagation  regimes  are  identified  and  analyzed:  a  fast  recombination 
regime  in  which  the  recombination  reaction  occurs  in  a  steady-state  manner  with  the  high-activation- 
energ\  branchtng  reaction  in  an  overall  reaction  region  which  is  much  thinner  than  that  of  diffusion,  and 
an  intermediate  recombination  regime  in  w  hich  the  recombination  region  is  much  thicker  than  the  branching 
region  but  much  thinner  than  the  diffusion  region.  The  analysis  yields  the  characteristic  dual  solution, 
e.xtinction  turning  point  flame  response,  and  shows  that  the  flame  propagation  speed  is  reduced  to  e' '  ' 
of  the  adiabatic  value  for  both  reaction  mechanisms.  The  generality  of  this  limit  flame  speed  is  noted. 


1.  INTRODUCTION 

A  CLAisic.xL  model  problem  for  the  study  of  premixed 
flame  extinction  is  that  of  Spalding  (!].  who  analyzed 
the  one-dtmcnsional  freely -propagating  flame  with 
temperature  sensitive  one-step  overall  reaction  and- 
radiative  heal  loss.  By  letting  these  processes  assume 
power-law  temperature  v  ariations,  with  the  exponents 
being  1 1  and  4  for  the  reaction  rate  and  heat  loss  rate 
such  that  the  former  is  more  temperature  sensitive 
than  the  latter,  the  analysis  yields  an  extinction  turn¬ 
ing  point,  at  which  the  normalized  flame  speed 
=  ii  f.  is  0.504.  Where  s,  is  the  flame  speed  and 
Sf  the  adiabatic,  laminar  flame  speed,  and  the  sub¬ 
script  'ex'  designates  the  extinction  state.  A  more  rig¬ 
orous  analysis  [2]  of  the  Spalding  problem  by  using 
activation  energy  asymptotics  and  for  weak .  0(c)  con¬ 
ductive  heat  loss  subsequently  yielded  j,.,,  =  e" '  % 
which  is  about  0.607.  and  a  suitably-scaled  heat  loss 
rate  L  «=  e‘ '  at  the  state  of  extinction,  where  t  is  the 
reciprocal  of  the  nondimensional  activation  energy  to 
be  defined  later.  The  same  result  was  shown  to  also 
hold  [3]  for  a  general  0{r.)  volumetric  heal  loss  func¬ 
tion. 

Recently.  Sibulkin  and  co-workers  [4-6]  reported 
numerical  solutions  of  transient  planar  and  out¬ 
wardly-propagating  lean  methane  air  flame,  allowing 
for  radiative  heat  loss  but  with  constant  transport 
properties  and  one-step  overall  reaction.  It  was  shown 
that  at  the  limit  of  propagation  i,  again  assumed  a 
Value  around  0.6.  Laksnmisha  el  al.  [7],  and  Law 


and  Egolfopoulos  [8]  extended  these  calculations  to 
include  variable  properties  and  detailed  chemistry, 
and  further  demonstrated  the  near-constancy  of  this 
value.  Compared  to  previous  analytical  studies,  these 
computational  results  are  significant  because  they 
indicate  the  possibility  that  at  the  extinction  limit  the 
flame  speed  is  always  reduced  to  about  60%  of  its 
adiabatic  value,  being  very  insensitive  both  to  the 
reaction  mechanism,  whether  it  is  one-step  or  detailed, 
and  to  the  loss  mechanism,  whether  it  is  conductive 
or  radiative,  as  long  as  it  is  0(e)  and  volumetric  in 
nature. 

The  objective  of  the  present  study  is  to  provide  a 
fairly  general  proof  of  the  above  possibility  based 
on  activation  energy  asymptotics.  Recognizing  that 
previous  analytical  studies  mostly  involve  a  one-step 
large  activation  energy  reaction,  which  obviously  can¬ 
not  capture  the  inherently-important  chain  branching 
and  termination  nature  of  realistic  reaction  schemes, 
we  shall  employ  the  Zerdovich-Linan  two-step  mech¬ 
anism  (9,  10],  which  consists  of  a  branching  reaction 
and  a  competing  recombination  (termination)  reac¬ 
tion.  This  is  believed  to  be  the  simplest  representation 
of  the  chain  nature  of  realistic  reaction  mechanisms. 
A  general,  0(e)  volumetric  heat  loss  function  will  be 
used  in  the  analysis. 

The  system  to  be  analyzed  is  formulated  in  the  next 
section.  It  will  be  shown  that  there  are  two  situations 
of  interest,  which  respectively  involve  fast  and  inter¬ 
mediate  recombination  reaction  rates.  These  two  situ¬ 
ations  are  separately  analyzed  in  Sections  3  and  4. 
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NOMENCLATURE 


a.  b.  c  integration  constants 

B  pre-exponential  factor 

Cp  specific  heat  at  constant  pressure 

D  mass  diffusion  coefficient 

Da  Damkdhier  number 

£  activation  energy 

/  mass  flux  through  the  flame  or  the  flame 

propagation  rate 

/  flame  propagation  rate  in  the  adiabatic 

limit 

F  rate  controlling  reactant  (fuel) 

H  temperature-dependent  volumetric 

heat  loss  function 
Heaviside  function  defined  after 
equation  (A24) 

L  heat  loss  coefficient 

Lc  conductive  heat  loss  coefficient 

£ii  radiative  heat  loss  coefficient 

Le  Lewis  number 

M  third  body  of  collision  needed  for 

recombination 
P  combustion  product 

p  pressure 

qv  heat  of  combustion  per  unit  mass  of  fuel 
R  radical  species 

R  universal  gas  constant 

5,  flame  speed 

S(  adiabatic,  laminar  flame  speed 

T  temperature 

adiabatic  flame  temperature 
T,  leading  order  flame  temperature 

H'  molecular  weight 

.V  spatial  coordinate  attached  to  the  flame 

front 

)■  mass  fraction. 

Greek  symbols 

r  parameter  defined  after  equation  (A20) 


V  Euler’s  constant 

d  small  parameter  representing  the 

relative  thickness  of  the 
recombination  region  to  the  preheat 
region 

e  small  parameter  of  expansion,  defined 

as 

tj  compressed  coordinate  defined  as 

£X 

©  temperature  perturbation  in  the 

recombination  region 
0  temperature  perturbation  in  the 

branching  region 
A  burning  rate  eigenvalue 

X  thermal  conductivity 

i  stretched  coordinate  defined  as  .f/e 

p  density 

d>  fuel  concentration  perturbation  in  the 

branching  region 

V  radical  concentration  in  the 
recombination  region 

^  radical  concentration  in  the  branching 

region 

(»  reaction  rate  function. 

Subscripts 

ex  the  extinction  state 

F  fuel 

M  third  body  of  collision  needed  for 

recombination 
R  radical  species 

—  3c  quantities  at  the  upstream  ambiance. 

Superscripts 

— .  -I-  quantities  in  the  upstream  and 

downstream  of  the  branching  region 
n  pressure  exponent 

~  dimensionless  quantities. 


2.  FORMULATION 

The  two-step  branching-recombination  reaction 
mechanism  adopted  in  this  study  is  given  by 

F-I-R-2R  (1) 

R  +  R+\t  P+  M  (2) 

where  F.  R.  P  and  M  respectively  denote  the  rate 
controlling  reactant  (say  the  fuel),  a  radical  species, 
the  combustion  product,  and  a  third  body  of  collision 
needed  for  recombination.  In  this  mechanism,  reac¬ 
tion  (I)  represents  the  two-body,  high  activation 
energy,  branching  reaction  that  has  a  small  heat  of 
combustion  and  produces  more  radicals  than  it  con¬ 
sumes  when  reacted  with  the  fuel.  Equation  (2)  is  the 
three-body,  low  activation  energy .  highly  exothermic, 
termination  reaction  that  combines  the  radicals  to 


form  combustion  products  and  generate  heat.  For 
simplicity  we  shall  therefore  assume  that  the  bran¬ 
ching  reaction  is  thermally  neutral  while  the  ter¬ 
mination  reaction  has  zero  activation  energy.  If  we 
further  assume  that  the  rates  of  these  two  reaction 
steps,  cij,  vary  with  the  first  order  of  each  reactant, 
then  they  can  be  respectively  expressed  for  reactions 
(1)  and  (2)  as 


w,  »  ^ exp  (-£//?  D, 


W;  =  BzP"- 


(3) 

(4) 


where  T  is  the  temperature.  >’,  and  H',  the  mass  frac¬ 
tion  and  molecular  weight  of  species ;.  p  the  pressure. 
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B,  and  n,  the  pre-exponential  factor  and  pressure 
exponent  of  reaction  /.  £  the  activation  energy  of  the 
branching  reaction,  and  R  the  universal  gas  constant. 
The  mass  fraction  of  the  third  body.  can  be  con¬ 
sidered  to  be  a  constant. 

The  physical  problem  under  stud\  is  the  steady 
propagation  of  the  one-dimensional  premixed  flame 
in  the  doubly-infinite  domain  with  an  arbitrary  tem¬ 
perature-dependent  volumetric  heat  loss  function 
H(T).  The  governing  equations  and  boundary  con 
ditions  are  given  by 

dr  d  /  dr\ 


.Y--x:  r-r, V, - 

„  d)V  ^ 

.V  -  X  ;  -*  0.  - - ►  0. 

d.v  d.Y 


)V  . , 

d.Y 


-0. 


.0. 


(8a) 

(8b) 


where  .y  is  the  spatial  coordinate  attached  to  the  flame 
front.  /  the  mass  flux  through  the  flame,  p  the  density. 
(Y  the  specific  heat  at  constant  pressure.  </^  the  heat  of 
combustion  per  unit  mass  of  fuel  consumed.  /.  the 
thermal  conductivity,  and  D,  the  mass  diffusion 
coefficient  of  species  /.  The  quantities  A.  and  pD, 
are  assumed  to  be  constants.  It  is  appropriate  to  con¬ 
centrate  on  only  0(r.)  heat  loss  because  it  represents 
the  limiting  situations  of  the  extinction  of  weakly- 
burning  flames.  Clearly  a  flame  will  extinguish  with. 
0(  1 )  heat  loss  if  it  extinguishes  w  ith  Oft)  heat  loss. 

It  is  worth  noting  that  equation  (8b)  implies 
T  -»  T , .  y,  -•  0.  and  )’r  —Oai  x  -»  x .  for  w  hich  T, 
is  the  adiabatic  flame  temperature  if  the  system  is 
adiabatic  and  T,  =  r_ ,  in  the  presence  of  heat  loss. 
Thus  both  of  these  conditions  are  applicable  and  yield 
the  same  result.  Equations  (8b)  are  adopted  because 
they  arc  more  general.  For  the  case  of  moderate  (i.e. 
0(1))  heat  loss  in  the  downstream  of  the  branching 
region,  equations  (8b)  are  still  valid  except  )>  ^  0  at 

.Y  -•  X  . 

Introducing  the  nondimensional  quantities 


f  = 


>% 


^  R 
£ 
Lc 


Yy  -r  »  V  ■  f  ■  i.(  fc,y 

2E  R  ^  _  2Hi 

.  ‘r'  (  f  )-tpV^  ' 


Du, 


B,y  y, 

( ()-( "ir,"  ■ 


iB,>.p'-  y,  .,  Fm 
(/  I'fpMVH’M  ■ 

where  Le,  is  the  Lewis  number  of  species  i.  Da,  the 
Damkohler  number  of  the  yth  reaction.  f  '(Lev)  the 
flame  propagation  rate  in  the  adiabatic  limit,  and 
f,a  =  1  -I-  f_,  the  adiabatic  flame  temperature,  the 
governing  equations  are  nondimensionalized  to 

(9) 


= -Da,Y^9^np(-£!T),  (10) 

J_^R  ^ 
d.v  Le„  d.Y- 

Z)o,rpfRexp(-£,f)-/)fl:fi.  (II) 
.Y-*-x:  f-f-,.  >V-I.  ?R-0.  ^-»0. 

(12a) 

df  „  df;  „  df’R  „  „ 

.Y  -»  X  :  trr  -*  0.  -  0.  -  0.  ^  -  0. 

d.Y  d.Y  d.Y 

(12b) 


This  system  w  ill  be  solved  by  activation  energy  asymp¬ 
totics. 

Depending  on  the  rate  of  the  recombination  reac¬ 
tion  relative  to  that  of  the  branching  reaction,  three 
flame  propagation  regimes  can  be  identified.  In  a  fast 
recombination  regime  both  reaction  rates  are  of  the 
same  order  so  that  the  reactions  occur  in  the  same 
thin  reaction  region.  In  an  intermediate  re¬ 
combination  regime  the  recombination  rate  is  much 
slower  than  the  branching  rate  but  is  much  faster  than 
the  diffusion  rate.  Consequently  the  recombination 
region  is  much  thicker  than  the  branching  region  but 
much  thinner  than  the  preheat  region.  Finally,  in  a 
slow  recombination  regime  the  recombination  rate  is 
either  comparable  to  or  slower  than  the  diffusion  rate 
such  that  the  thickness  of  the  recombination  region  is 
either  of  the  same  order  of  or  larger  than  that  of  the 
preheat  region. 

In  the  next  two  sections,  the  fast  and  intermediate 
recombination  regimes  will  be  analyzed  sequentially. 
The  slow  recombination  regime  will  not  be  analyzed 
because  it  represents  very  weak  chemical  systems 
w  hich  are  not  likely  to  be  of  interest  to  combustion. 


3.  FAST  RECOMBINATION  REGIME 

In  this  regime,  the  rates  of  the  recombination  and 
branching  reactions  are  of  the  same  order  so  that  both 
reactions  occur  in  the  same  thin  reaction  region.  Thus 
the  radicals  are  consumed  to  produce  heat  almost 
immediately  after  they  are  generated.  This  implies  that 
the  concentration  of  the  radicals  is  basically  inde¬ 
pendent  of  time  and  is  very  small,  and  hence  the  steady 
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siaieapproximation.d)  ,,  d/  -►  0.  is  applicable.  Figure 
1  shows  the  species  and  temperature  profiles  of  such 
a  flame. 

With  the  application  of  the  steady  state  approxi¬ 
mation.  setting  equation  ( 1 1 )  to  zero  readily  yields 

?R  =  ^fvexp(-£f).  (13) 

The  vanishing  of  the  LHS  of  equation  (11)  is  justified 
by  the  small  radical  concentration.  Substituting  Pr 
into  the  reaction  rate  terms  of  equations  (9)  and  (10). 
the  problem  is  then  reduced  to  a  one-step  reaction 
with  a  rate 

gf^exp,-2^f). 

w  hich  is  second  order  with  respect  to  .  and  has  an 
effective  Damkdhier  number  Da]  Da^  and  an  effective 
activation  energy  2£.  Compared  with  previous 
asymptotic  studies  (see.  for  example,  refs.  (2, 3])  which 
adopted  a  one-step  first-order  overall  reaction,  the 
difference  in  the  mathematical  aspect  of  the  two  prob¬ 
lems  is  only  the  reaction  order.  It  is  therefore  of  inter¬ 
est  to  study  that,  apart  from  the  obvious  difference  in 
the  physical  interpretation  of  these  two  problems, 
in  terms  of  the  chain  branching  and  termination 
reactions  through  modifications  of  the  effective 
Damkdhier  number  and  activation  energy,  what 
additional  effect  a  simple  change  in  the  reaction  order 
can  produce.  Due  to  the  similar  nature  of  the  asymp¬ 
totic  derivation  with  previous  studies,  only  the  key 
steps  w  ill  be  show  n  below . 

In  the  outer,  chemicalh  inert  regions,  there  is  no 
branching  reaction  in  the  upstream  region  because  the 
flow  temperature  is  low  and  the  reaction  has  a  high 
activation  energy.  In  the  downstream  region  branch¬ 
ing  is  terminated  because  of  complete  fuel  consump¬ 
tion.  The  recombination  reaction  does  not  exist  in  the 
outer  regions  due  to  complete  radical  consumption. 
Thus  the  outer  solutions  for  fv  are 

=  1 -(h„-l-£h, -f-  ■)cxp(LCf/v).  (14) 
)'V=0.  (15) 


Fio  I.  Scheniuiic  of  the  flame  siruciure  for  the  fast  reconr 
binjlKin  regime,  has  been  magnified  forclarin. 


while  f  can  be  expanded  as  f-  =  fg  H — 

This  yields 


fi  =  ^-x  -l-floCxp(/v), 

(16a) 

(16b) 

(17) 

In  the  above  the  superscripts  —  and  -f  respectively 
denote  quantities  in  the  upstream  preheat  region  and 
downstream  burned  region,  while  a,  and  b,  the  inte¬ 
gration  constants  to  be  determined.  Note  that  in  the 
downstream  region,  a  compressed  coordinate  »  ex 
is  used  because  the  temperature  decreases  gradually 
from  the  flame  temperature  to  the  ambient  tem¬ 
perature  due  to  the  heat  loss. 

In  the  inner,  chemically  reactive  region,  the  co¬ 
ordinate  is  stretched  as  {  =  x/c  while  f  and  are 
expanded  as 


f,  —  c6,—c-62  +  --. 

(18) 

=  £^1,  - 

(19) 

Substituting  equations  (18)  and  (19)  into  equations 
(9).  ( 10)  and  ( 1 3),  and  expanding,  we  obtain  the  inner 
equations 

=  A.^,exp|^ 

(20) 

.  1  d=«, 

d;'  Le^  di= 

(21) 

d-’O;  1  d->:  Jde,  ddiA 

d;-  U,  ^\di  di)~ 

(22) 

where 

(23) 

Matching  the  inner  and  outer  solutions  yields  the 
relations 

Oil  =  fi-  f-,  ,  fn  (»1  =  0)  =  ff,  bo  « 

1,  (24) 

and  the  matching  conditions  to  solve  equations  (20) 
to  (22). 

6,(i-  -x)=  -fr(.v  =  0)-/(fr-f.x);. 

(25a) 

0,(^— «)  =  -fr(»/  =  o). 

(25b) 

-  fHT-o)dx, 

(26a) 

df<.  -  ,  df:  /?(ff) 

(26b) 

<t>,U-  -x)  >=  -b,-LeJi. 

(27a) 
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d>i(:-y.)  =  0. 

(27b) 

(28a) 

d<f>- 

^(i-^-)  =  0. 
d^ 

(28b) 

where  equations  (27a.  b)  also  imply 

^  (i  -  y )  =  -  Le^f. 

(29a) 

(s-«y)  =  0. 
d. 

(29b) 

Equations  (16b)  and  (17)  have  been  applied  in  obtain¬ 
ing  the  second  relation  of  equations  (26a.  b). 

Integrating  equations  (21 )  twice  and  (22)  once  sub¬ 
ject  to  the  matching  conditions  n  equations  (2S)-(28) 
yields 

f,  =  1+f.,  = 

(30) 

(31) 

where 

£  =  2|^/?(f.j+/j'  /?(/,; )d.vj. 

(32) 

It  can  be  readily  demonstrated,  by  re-defining 
.V  =  fx  =  .V  [/.  (/c>)].  that  £  is  independent  of  / 
Equation  (30)  shows  that  the  leading  order  flame 
temperature  is  the  adiabatic  flame  temperature,  which 

is  reasonable  because  the  heat  loss  is  assumed  to  be 
C7(e).  Substituting  equation  (31)  into  equation  (20) 
and  integrating  the  resulting  equation  once  subject  to 
the  boundary  conditions  in  equations  (27a.  b)  and 
(29a.  b).  we  obtain  an  expression  which  determines 
the  flame  propagation  rate. 

7=exp(^^)  =  ^A.  (33) 

In  the  adiabatic  limit.  /?  =  0  and  /  =  I .  We  then  have 
(i.e^A  2)  «  I  for  the  laminar  flame  propagation  rate 
/  .  Equation  (33)  then  becomes 


The  quantities  f  f  and  />,  are  not  of  interest  and  hence 
w  ill  not  be  solsed. 

Equation  (.34).  of  course,  is  exactly  the  same  as 
that  of.  sa>.  Joulin  and  Clavin  [2]  who  considered 
conductive  heat  loss  and  a  hrst-ordcr  reaction.  Repro¬ 


Fic.  2.  Normalixed  flame  propagation  rate  /as  a  function 
of  the  heat  loss  parameter  L  for  the  fast  and  intermediate 
recombination  regimes. 

ducing  equation  (34)  in  Fig.  2,  it  is  seen  that  there 
exists  a  maximum  value  of  £.  above  which  no 
solution  exists.  For  £  <  £„,  there  are  two  solutions 
for  each  £.  although  it  is  well  established  that  only 
the  upper  branch  gives  the  stable  solution.  This  critical 
state  is  then  defined  as  that  of  extinction,  being  charac¬ 
terized  by 

£,,  =  e- '  and  /„  =  e* '  /  (35) 

4.  INTERMEDIATE  RECOMBINATION  REGIME 

In  this  regime  the  rate  of  the  recombination  reaction 
is  much  slower  than  that  of  the  branching  reaction 
but  much  faster  than  the  diffusive-convective  trans¬ 
port  rate.  Consequently  the  recombination  region  is 
much  thicker  than  the  branching  region  but  much 
thinner  than  the  transport  region.  By  defining  a 
second  small  parameter  6  to  describe  the  characteristic 
thickness  of  the  recombination  region,  where  e  « 
1,  the  flame  structure  then  includes  an  0(e) 
branching  region  sandwiched  by  0{S)  recombination 
regions,  which  in  turn  are  embedded  within  an  0(1) 
upstream  preheat  region  and  an  0(1 /e)  downstream 
burned  region.  Since  the  exothermic  recombination 
reaction  continues  subsequent  to  completion  of  the 
branching  reaction,  temperature  will  also  con¬ 
tinuously  increase  until  heat  loss  becomes  dominant. 
The  species  and  temperature  profiles  are  shown  in 
Fig.  3. 

B\  separately  analyzing  the  five  regions  and  per¬ 
forming  the  requisite  matching,  which  is  presented  in 
the  Appendix,  the  flame  speed  response  is  now  given 
by 


Comparing  equations  (36)  and  (34).  we  can  see  that 
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Fig.  3.  Schematic  or  the  flame  structure  for  the  intermediate 
recombination  regime. 


their  functional  expressions  are  similar  although  the 
branching  and  recombination  reactions  occur  at 
different  length  scales.  Thus  the  behavior  of  /  as  a 
function  of  L  is  qualitatively  similar  to  that  of  the  fast 
recombination  regime,  as  shown  in  Fig.  2.  The  flame 
extinction  limit  is  now  described  by 

7„=e-'=.  4.  =  (4  3)c-‘.  (37) 

It  is  interesting  to  note  that  although  the  reaction 
kinetics  adopted  in  this  regime  are  different  from  those 
of  Section  3.  the  flame  still  extinguishes  at  /=  e" '  % 
albeit  at  a  different  L,^.  Combining  this  result  with 
those  obtained  from  the  analysis  by  adopting  a  one- 
step  overall  reaction  (2)  and  the  numerical  studies  {4- 
8].  it  may  be  suggested  that  extinction  is  achieved 
when  the  inherent  nonadiabaticity  of  the  system 
reduces  the  flame  propagation  rate  to  about  e' ' '  or 
0.607  of  its  adiabatic  value.  The  fact  that  a  large  heat 
loss  is  needed  to  extinguish  the  flame  in  the  inter¬ 
mediate  recombination  regime  as  compared  to  that 
in  the  fast  recombination  regime  is  also  reasonable 
because  the  overall  reaction  rate  is  less  temperature 
sensitive  for  the  slower  recombination  rate. 

5.  CONCLUDING  REMARKS 

In  this  stud) .  we  have  analyzed  the  extinction  limit 
of  the  planar  premixed  flame  with  0(e)  volumetric 
heat  loss  by  adopting  a  two-step,  branching  and  ter¬ 
mination  kinetics.  Based  on  the  relative  rates  of  these 
two  reactions,  the  fast  and  intermediate  combination 
situations  are  studied.  The  results  show  that  the  flame 
always  extinguishes  t^hen  the  flame  propagation  rate 
is  reduced  to  e" ' '  of  its  adiabatic  value,  which  is 
consistent  with  previous  analytical  results  with  one- 
step  reaction  and  numerical  results  with  detailed 
chemistry. 

The  abo\e  analysis  is  based  on  a  general  heat  loss 
function  H.  Two  loss  mechanisms  that  are  usually 
specified  are  conduction  and  radiation.  For  a  linear 
heat  loss  function.  H  «  Lr(  T-  T. , )  where  4  is  the 
heat  loss  coefficient.  3\e  then  have  L  =  (A/.Lc) 
i  f  f^)^  For  radiative  loss.  H  Lk(T* -  T*. ,  )  and 
we  have 


Finally,  we  note  that  the  influence  of  the  reactant 
concentration  .?n  flame  extinction  is  primarily 
through  the  factor  (/')*  in  the  definition  of  fi.  Thus 
as  the  fuel  concentration  becomes  either  leaner  or 
richer,  the  decrease  in  p  would  lead  to  a  cor¬ 
responding  increase  in  R.  It  is  therefore  clear  that  the 
present  result  also  predicts  the  existence  of  con¬ 
centration  limits  at  the  heat  loss  rate  £„,  beyond 
which  steady  flame  propagation  is  not  possible.  Such 
limits  have  been  identified  as  the  flammability  limits 
[1,4,5.  7,8]. 
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APPENDIX :  DERIVATION  FOR  THE 
INTERMEDIATE  RECOMBINATION  REGIME 

All  the  solutions  are  expanded  with  respect  to  the  two 
small  parameters  c  and  6.  Thus  in  the  outer,  transport  regions 
away  from  the  recombination  region,  ffi  »  0.  while  Tf  is 
given  by  equation  (15)  and 

ff  =  1  —  I[/>i,-t-0(6)]  +  £l/)| +  0(^1) 

-»-0(£-);  expfLfr.^^)  (AI) 
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The  solution  of  T-  can  be  obtained  b\  first  expanding  it  in 
the  form 

T-  =  Iff  +Sfi  +0(6=)J+tlf5  +0((5)]-^0(c=).  (A2) 

which  is  then  substituted  into  equation  (9)  to  >'ield  equations 
(16a).  (17).  and 


H-f  -4(r;+ff)=. 
4-?  =rc5/(r;+ff)’. 


f;  =a,exp(7i). 


>?(fo-)d.v. 


The  compressed  coordinate  i)  =  c.v  is  again  used  in  the 
burned  region  such  that  all  the  terms  in  f  *  are  functions  of  q. 

In  the  0{S)  recombination  region,  the  coordinate  is  stret¬ 
ched  as  ;  »  .V  6.  There  is  no  branching  reaction  in  this  region 
because  of  its  high  activation  energy.  Thus  the  outer  solution 
of  ff  is  still  applicable.  Moreover,  only  0(d)  vanations  are 
possible  for  f  -  and  f  g  so  that  their  expansions  are 

f-  =  -S-Qj  -t-0(<l-)] 

-cl©.'  -*-<5©i  -t-Old-tJ-t-Otc’).  (A5) 
f:  =  |d4<=  4.d=4<j  +0(d')] 

H-cCff -t-df; -t-0(^-)J-H0(c-).  (A6) 

Substituting  equations  ( A5)  and  ( A6)  into  equations  (9;  and 
(II).  but  with  the  branching  reaction  frozen,  and  expanding, 
we  obtain  for  f* . 

I  d-f^ 

- rrf- =  2'-'A;4'f4'^  1=1.;.  (A7) 


I  d-'4'=:  .df^ 
Lf*  d,'-  '  d; 


=  :A;I4'F4',1; 


-t-((-l)4'rH',--i].  1=1.2.  (A8) 

where  A;  =  S^Dn;  is  the  reduced  Damkdhler  number  for  the 
termination  reaction.  The  local  coupling  functions  are  given 

b> 

The  volumetric  heat  loss  is  not  important  in  the  above  equa¬ 
tions  because  the  recombination  regions  are  still  very  thin, 
of  Old).  Each  of  the  above  equations  needs  to  be  solved 
separately  in  the  upstream  and  downstream  of  the  branching 
region. 

Matching  the  solutions  in  the  recombination  region  with 
the  outer  diffusive-convective  regions  yield  a,.  =  f,  o  —  7.  , . 

and  the  conditions  to  solve  equations 

(A7)-(A10). 


4’,-  (,'  -  ±  X )  -  0.  / 
e;  (,'  -.  T. )  =  -  T:  (»;  - 
d0i ,  df,: , 


1.2. 3. 4. 


©7(F-»-x)«- 
©F  (J—  -  X) 


•ti  —  7. , ),. 

-  -  f :  (.V  -  0). 


-j*  fl(f,:)dv.  (AI6) 

Equations  (17)  and  ( A4 )  have  respeclivelv  been  used  in  deriv¬ 
ing  the  second  relation  of  equations  (.AI3l  and  <AI6). 

Solving  equations  (.A7i  and  (.A8)  subject  to  the  matching 
conditions  in  equation  i  A 1 1 ).  we  obtain 


rc-t-cf  [(rc-(-cf)-'"5  r^]’ 

4^  r  r  fi  3  cfcfr  1  J 
d;  '(rc-hffi’Lrc+cr  2(r;-fcf)-  s"- 


where  c,-  are  the  integration  constants  to  be  determined  and 
r  =  (2Le*A:'3)'=. 

Applying  equations  (AI7)-(A20)  to  equations  (A9)  and 
(AIO).  then  solving  the  resulting  equations  subject  to  the 
matching  conditions  (A12)-(AI6).  we  obtain 


-(7(fr.o-f-x)C-n.)«.(-;)-fr(';  =  0)H.(;).  (A21) 
©-  »  Tf? 

'•  ifR  (r;-i-rf )’ 

-  7r  (.V  =  0)/f.  ( -  0  -  7?  (n  =  0)H,  C).  (A22) 
d=©j  8  [■  ctA:  ,/6  1  \'] 

d:-  (r;-hcf)’L(r;+cf)='^'-'l.5  ujj 

-7'(7,.o-f-,)H.(-C).  (A23) 

d©i  1  r  r  rj  ,V6 

d;  (TC-^cf )’  rc -f-cr  ~  Vj ' 

fyff  A;  1  ^(ft.o)  ,, 

-(TTTTfpJ'-r"''-’ 

-r/fc(,:-o)-r  (A24) 

where  H,  is  the  Heaviside  function  with  3/,(C)  *  1  for  J  >  0 
and  0  for  V  <  0.  In  the  above  only  the  solutions  required  for 
matching  with  those  in  the  branching  region  are  presented. 
Because  (d4'r  id;)o  >  0  and  (d4'F/dC)o  «  0.  we  must  have 
cF  <  Oand  rf  >  0, 

In  the  thin  branching  region,  the  coordinate  is  stretched 
as  i  =  .x/t  while  ff  and  f*  are  expanded  as 

ff  =  £(<f)|-i-  ••)-(-r-(^;-l-  '  )-l-0(£’).  (A25) 

=  {Sill,  +  S' ill, -i - )-g£(l/'.--l - )-(-£'(|f>j-t--  •  •)-(-0(£’). 

(A26) 

Because  the  activation  energy  for  the  branching  reaction  is 
high,  variation  of  temperature  can  only  be  0(e).  Thus  f  is 
expanded  as 

f  =  (7(0  — 67t.i  -t — )  — c’(®3  +  '  ■  Oh - 


Although  the  activation  energy  is  large,  it  is  not  extremely 
high  such  that  the  branching  reaction  can  occur  in  a  tem¬ 
perature  range  smaller  than  its  maximum  value  by  an  0(6) 
amount.  Substituting  equations  (A25)-(A27)  to  equations 
(9)-(II)  and  expanding,  we  obtain  the  following  equations 
that  describe  the  branching  region. 

^=0.  1=1.3.  (A28) 


A  /  ' 
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6-0.  ydO, 

-,-TV  -T  -TT  = 

d<-  d,- 

d--\Lcy  UtJ 


(A31) 

(A32) 


(A45) 


1  d-0;  ,d(<>,  1  d<lly  y.dV: 

—  — ^ 4- —  — ^ IA33) 
Ley  d,"  d^  Lct  a;-  a, 

1  d'di 


where 


M  d.:--  (A34) 

A,=t^dO^,exp(-y,--^|-— )  (A35) 

is  the  reduced  Damkdhier  nurnber  for  the  branching  reac¬ 
tion. 

Matching  is  performed  first  b\  expressing  the  outer  solu¬ 
tions  of  in  terms  of.;  =  .f  c.  expanding,  and  then  equating 
the  resulting  expression  with  the  solution  in  the  branching 
region.  This  gives  fco  *  1  and  equations  (27)-{29).  Next,  the 
solutions  of  f  =  and  fg  in  the  recombination  region  are 
expressed  in  terms  of  =  (c  S)~  and  expanded  to  yield 


Substituting  equations  (A44)  into  equation  (A36),  Ty,  can 
be  determined  if  f  J  (i;  =  0)  is  known.  Because  Ty,  i  represents 
the  temperature  at  which  branching  reaction  occurs,  it  is 
only  a  function  of  activation  energy  and  chemical  reactivity 
but  does  not  depend  on  the  heat  loss.  Thus  by  knowing  that 
(»}  =  0)  =  0  and  /=  I  in  the  adiabatic  limit,  we  obtain 


(A46) 


and  the  matching  conditions 


p-fr(n  =  0). 

(A-76) 

tion  (A34).  and  defining  thi 

r 

4>  =  T-, — . 
lU, 

:  1  =  4  (« f )’. 

(A37) 

Hcj  -8cl  (c;  )’. 

f  9^ 

<A38) 

r nW 

r  f  i 

lt<f)--'"5  r  J- 

(A39) 

and  hence  (i;  *  0)  =  fi.tif'  '  —  1). 

Next,  by  sequentially  integrating  equations  (A28)  and 
(A39)  twice,  and  equations  (A30).  (A3I)  and  (A33)  once, 
subject  to  equation  (A44)  and  the  proper  boundary  con¬ 
ditions  of  (A39)-(A43).  we  obtain  J'f.o®  l  +  f-« 

=  Cj .  r;  =  cj ,  hi  =  -  Ley  Jet  /8.  as  well  as 

where  L  is  defined  in  equation  (32). 

Finally,  substituting  equations  (A4S)  and  (A47)  into  equa- 


7= ’exp 


Gt)] 


(A48) 


dti4,.  ,  r  (ci 


±x)  = 


;  -8;)-{f:(.\  =0) 


LcrIcF I 

+J(fu-  f. ,  )cl/f.<-<>-  7r<';  =  0)/f.(c).  (A4I) 
O'M  -  i  y.)  =  2A;(4.;-<  ;-  ),=  (<-;)*.  <A42) 


dO,  . 

aT<v-'±x’-l  , 

d;  <c 


I  l<ir  [rjA: 

f)’|<rLf,  Ltrrl’ 

-j^  /■•■(  f,v  -  T. ,  >■  ^Jf:  ( f  =  0)  -I- 1  i9(  f o' )  d,vj«,  ( - 

(A43) 

In  analyzing  the  branching  region,  first  equations  (A28) 
for  \li ,  and  (A32)  are  solved  subject  to  equations  (27).  (A37) 
and  (A38)  to  yield 

(8A 

rLcV?:)  • 


(he  structure  equation  is  obtained. 

-d<J  -- - 

(A49) 

dv'  di, 

=  0)  =  In  ||^2LeFA,^^^^  J 

X  1^7'’ exp (^)]  j-21n(f'2).  (A50) 

<^(i-vx)  =  0.  (A5I) 

This  is  a  rnodified  Bessel  function  of  order  zero  whose 
solution  is  ^  s  cKM  )•  where  c  is  the  integration  qpnsunt, 
after  ^uation  (A51 )  is  applied.  Because  KM  »  0)  — 
-  In  (i  2)  -  )•.  where  •/  =  0.5772  is  the  Euler's  constant,  upon 
applying  equation  (ASO)  we  obtain 


|^2L(>fA,|^^^  /‘’exp^^jj  »exp(-2)). 


(A52) 


In  the  adiabatic  limit.  7  =  0  and  7=  I  •  we  have 

[^2LcfA,(^J  ’j  -  exp(-2v).  (A53) 

Substituting  equation  (AS3)  in  equation  (A52),  we  obtain 
the  flame  response  given  by  equation  (36). 
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The  Structure  of  Premixed  Methane-Air 
Flames  with  Large  Activation  Energy 
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We  examine  the  structure  of  both  lean  and  stoichiometric  premixed  methane-air  flames  with  a  reduced 
reaction  mechanism.  Our  starting  point  is  a  four-step,  C|-chain  mechanism  that  has  been  derived  previously 
using  a  series  of  steady-state  and  partial  equilibrium  assumptions.  This  same  mechanism  has  been  adopted  in 
several  recent  asymptotic  studies  that  have  used  the  ratio  of  the  branching  to  propagating  reactions  as  the 
perturbation  parameter  to  analyze  the  fuel  consumption  zone  within  the  flame  structure.  In  the  present 
study,  we  assume  that  the  aaivation  energS'  of  the  intermediate  reactions  in  the  fuel  consumption  zone  are 
sufficiently  large  to  employ  the  method  of  large  activation  energy  asymptotics.  We  obtain  temperature  and 
species  profiles,  as  well  as  a  structure  problem  whose  solution  determines  the  burning  rate  eigenvalue  in 
terms  of  a  parameter  that  represents  the  ratio  of  branching  to  terminating  reactions.  When  this  parameter  is 
set  to  zero,  the  results  of  previous  rate-ratio  analyses  are  recovered.  In  the  opposite  limit  that  this  parameter 
becomes  large,  the  structure  reduces  to  that  of  Linin's  premixed  burning  regime.  In  both  limiting  cases,  we 
determine  the  parametric  dependence  of  the  burning  rate  on  equivalence  ratio,  pressure,  and  the  ratio  of 
competing  rates  of  branching  and  termination  reactions.  Hie  trends  are  found  to  be  largely  in  agreement 
with  experimental  obsersations.  One  advantage  of  the  present  approach  is  that  the  exponentially  nonlinear 
reaction  rate  terms  are  retained,  thus  permitting  the  study  of  the  response  of  methane-air  flames  to  small 
perturbations. 


NOMENCLATURE 

P 

pressure 

C* 

non-dimensional  heat  re¬ 

A 

parameter  defined  in  Eq.  4.9 

lease  of  global  reaction  k, 

Aj,  Ai, 

scaled  preexponential  fac¬ 

defined  in  Eq.  2.4 

tors  defined  after  Eq.  4.5 

.  9 

parameter  defined  in  Eq.  S.5 

Aj,  Ajj 

preexponential  factors  in  re¬ 

9k 

heat  release  per  mole  of  fuel 

action  1,  11 

consumed 

integration  constants 

9f-9h-9h,->9co 

parameters  defined  after  Eq. 

integration  constants 

4.10 

‘■p 

specific  heat 

/?" 

universal  gas  constant 

D 

Damkohler  number  defined 

f 

temperature 

in  Eq.  4.10 

t 

temperature  perturbation  in 

D 

scaled  Damkohler  number 

oxidation  layer  defined  in 

defined  in  Eq.  6.21 

Eq.  5.2 

D, 

mass  diffusivit}'  of  species  i 

V 

velocity 

E^,  E^^ 

activation  energy  of  reac¬ 

W 

variable  defined  in  Eq.  4.18 

tions  1,  11 

molecular  weight  of  species 

F 

fuel 

1 

K 

equilibrium  constants 

overall  reactions  defined  in 

k. 

elementary  rate  constants 

Eq.  2.1 

Lc, 

Lewis  number  of  species  i 

X, 

sp-icies  variable  related  to 

m 

parameter  defined  in  Eq. 

■ 

mole  fraction  normalized  by 

(6.21) 

mole  fraction  of  fuel 

X 

spatial  variable 

Y 

scaled  variable  defined  in  Eq. 

*  Corresponding  author. 
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I', 

mass  fraction  of  species  i 

y. 

mole  fraction  in  fuel  con¬ 
sumption  layer,  defined  in 
Eq.  4.4 

z. 

mole  fraction  in  oxidation 
layer,  defined  in  Eqs.  5.2, 5.3 

Greek  Symbols 

a 

parameter  defined  in  Eq.  2.3 

y 

parameter  defined  in  Eq.  6.2 

8 

thickness  of  fuel  consump¬ 
tion  layer,  defined  in  Eq.  4.3 

e 

thickness  of  oxidation  layer, 
defined  in  Eqs.  5.10,  5.21 

{ 

stretched  spatial  coordinate 
in  fuel  consumption  layer 

V 

stretched  spatial  coordinate 
in  oxidation  layer 

V, 

Chaperon  efficiency  of 
species  i 

e 

temperature  perturbation  in 
fuel  consumption  layer 

A 

thermal  conductivity 

M 

thickness  of  radical  con¬ 
sumption  layer,  defined  in 
Eq.  4.21 

scaled  spatial  coordinate  de¬ 
fined  in  Eq.  6.14 

p 

densit)’ 

a 

stretched  spatial  coordinate 
in  radical  consumption  layer 

T 

nondimensionai  tempera¬ 
ture 

equivalence  ratio 

nondimensionai  overall  re¬ 
action  rates 

Subscripts  and  Superscripts 

b 

burned  value 

1 

species.  /  »= 

f,C0,C0,,H,H,,H20,02 

k 

global  reaction  step 

0 

value  at  radical  consumption 
zone 

u 

unbumed  value 

1.  INTRODUCTION 

The  theorj'  of  flame  propagation  has  been 
greatly  advanced  during  the  past  20  years 
through  the  use  of  large  activation  energy 


asymptotics  to  simplify  the  nonlinear  reaction 
rate  terms  in  chemically  reacting  flows.  By 
modeling  the  complex  chemical  kinetics  to 
consist  of  a  single  overall  reaction  with  large 
activation  energy  confined  to  a  narrow  region 
in  the  flow  field,  a  mathematical  treatment  of 
the  governing  equations  is  often  possible. 

While  the  one-step  chemistry  model  facili¬ 
tates  mathematical  analysis,  it  is  recognized 
that  this  simplification  is  incapable  of  provid¬ 
ing  a  satisfactory  explanation  for  phenomena 
that  depend  largely  on  reaction  intermediaries 
and  thus  an  accurate  description  requires  a 
more  detailed  chemical  kinetic  scheme.  Since 
full  kinetic  mechanisms  are  frequently  too  large 
to  permit  either  mathematical  or  numerical 
analysis,  there  has  been  a  great  deal  of  recent 
work  devoted  to  identifying  the  most  important 
reactions  for  various  fuel -air  mixtures  in  order 
to  reduce  the  number  of  steps  to  a  point  where 
the  governing  equations  become  amenable  to 
analysis.  For  example.  Miller,  et  al.  [1]  have 
derived  a  “short”  mechanism  for  methane-air 
flames  consisting  of  about  40  steps  using  only 
C, -hydrocarbon  species.  The  resulting  ^stem 
of  equations  has  been  solved  numerically  to 
compute  burning  velocities  for  lean  and  stoi¬ 
chiometric  mixtures  over  a  wide  range  of  con¬ 
ditions  [2]. 

Although  these  “short”  mechanisms  repre¬ 
sent  a  significant  simplification  over  the  origi¬ 
nal  full  system,  the  number  of  equations  is  still 
much  too  large  for  a  mathematical  treatment. 
Furthermore,  in  analyses  based  solely  on  nu¬ 
merics,  it  is  frequently  difficult  to  identify  the 
most  important  parameters  that  mostly  influ¬ 
ence  the  global  properties.  Thus,  in  order  to 
bridge  the  gap  between  the  one-step  kinetics 
and  the  large  multistep  schemes,  Peters  [3]  has 
systematically  reduced  the  existing  “short” 
mechanism  down  to  a  mechanism  consisting  of 
four  overall  reactions  using  a  series  of  steady- 
state  and  partial  equilibrium  assumptions.  The 
resulting  four-step  model  involves  just  seven  of 
the  original  species  and  includes  kinetic  data 
from  only  five  elementary  reactions.  On  the 
basis  of  this  reduced  mechanism  the  basic 
structure  of  methane -air  flames  has  been  stud¬ 
ied  asymptotically  [4-7]. 

In  many  of  these  studies,  the  activation  en¬ 
ergies  of  the  individual  elementary  reactions 
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are  not  considered  to  be  large  enough  to  em¬ 
ploy  large  activation  energy  asymptotics.  The 
reaction  rates  are  assumed  to  have  a  power  law 
dependence  on  temperature,  rather  than  the 
conventional  Arrhenius  form,  and  the  ratios  of 
the  individual  reaction  rates  are  used  as  small 
parameters  in  the  asymptotic  analysis.  This  type 
of  approach  has  been  used  in  two-,  three-,  and 
four-step  models  to  calculate  burning  velocities 
over  a  wide  range  of  equivalence  ratio,  pres¬ 
sure,  and  temperature.  However,  the  power-law 
dependence  on  temperature  does  not  exhibit 
the  same  kind  of  sensitivity  to  temperature 
variations  as  the  Anhenius  approximation. 
Cbnsequently,  the  response  of  the  flame  to 
various  external  perturbations  is  not  veiy  sig¬ 
nificant. 

In  contrast  to  previous  works,  the  key  termi¬ 
nation  reaction  in  these  studies  was  found  to 
be  the  propagating  step,  rather  than  the 
three-body  recombination  reaction.  Further¬ 
more,  the  leading  order  flame  structure  analy¬ 
sis  yielded  a  value  for  the  “Damkohler  num¬ 
ber”,  which  was  independent  of  the  burning 
rate.  The  resulting  expression  was  therefore 
used  to  calculate  the  temperature  f°  at  the 
leading  edge  of  the  flame  zone.  A  solution  for 
f°  was  previously  sought  of  the  form  f  ®  ~  fj, 
-t-  0(c),  where  Tj,  is  the  known  flame  temper¬ 
ature  and  e  (assumed  small)  is  a  measure  of 
the  flame  width,  which  in  turn  is  proportional 
to  the  burning  rate.  Thus  an  expression  for  the 
burning  rate  was  obtained  by  equating  the  first 
two  terms  in  this  asymptotic  expansion  and 
solving  for  c. 

In  the  present  study,  we  adopt  the  four-step 
mechanism  as  a  starting  {>oint  to  analyze  the 
structure  of  premixed  methane -air  flames  in 
the  same  manner  as  Peters  and  Williams  [4]. 


That  is,  we  assume  that  the  H  radical  is  in 
steady  state  and  that,  to  first  approximation, 
the  water-gas  shift  is  in  partial  equilibrium, 
effectively  reducing  the  four-step  mechanism 
to  a  two-step  mechanism.  However,  as  in  [4], 
we  shall  incorporate  nonequilibrium  effects  as 
a  perturbation  from  the  leading  order  partial 
equilibrium  state.  Our  analysis  differs  from 
theirs  in  that  we  retain  the  Arrhenius  form  for 
the  individual  reactions  and,  except  for  the 
three-body  recombination  reaction  with  zero 
activation  energy,  we  consider  the  effective 
activation  temperatures  to  be  sufficiently  large 
such  that  activation  energy  asymptotics  can 
provide  reasonable  and  useful  results.  We  re¬ 
mark  that  although  activation  temperatures 
have  been  estimated  to  be  on  the  order  of  7(X)0 
K,  it  is  known  that  even  modest  values  are 
sufficient  to  confine  the  reaction  to  narrow 
regions,  thereby  enabling  the  use  of  activation 
energy  asymptotics  to  obtain  remarkably  accu¬ 
rate  results  [8]. 

In  our  study  the  flame  structure  remains 
identical  to  that  analyzed  by  Peters  and 
Williams  [4]  (Fig.  1).  It  consists  of  a  chemically 
inert  preheat  zone  where  the  temperature  in¬ 
creases  from  its  unbumed  value  to  a  character¬ 
istic  temperature  at  which  reaction  can  take 
place.  There  is  a  very  thin  fuel  consumption 
zone  where  the  hydrogen  reaction  occurs,  and 
the  H  radical  is  entirely  consumed  by  the  fuel 
in  a  yet  thinner  zone  that  defines  the  leading 
edge  of  the  flame  structure.  Immediately 
downstream  of  the  fuel  consumption  layer  lies 
a  relatively  broad,  but  still  asymptotically  small, 
oxidation  layer  where  Hj  and  CO  oxidize  to 
form  HjO  and  CO2.  This  entire  structure  is 
embedded  in  a  flow  field  which  is  chemically 
inert  due  to  the  absence  of  radicals. 


Fuel 

Consumption 
I— »  —I 


OsiOetion  Layer 


Fig.  1.  Schematic  of  the  flame  structure,  illustrat¬ 
ing  several  species  profiles  and  the  temperature 
distribution. 
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The  present  analysis  provides  a  complete 
description  of  the  above  flame  structure,  and 
the  role  of  each  of  the  overall  reaction  steps  is 
determined.  Our  results  include  the  evaluation 
of  the  burning  rate  in  terms  of  pressure,  equiv¬ 
alence  ratio,  and  a  parameter  that  represents 
the  ratio  of  the  competing  rates  of  the  branch¬ 
ing  and  termination  reactions.  In  particular,  we 
find  that  when  this  rate-ratio  is  reduced,  the 
burning  rate  eventually  decreases  to  zero.  For 
the  opposite  limit  that  the  rate-ratio  is  large, 
we  find  that  the  burning  rate  asymptotes  to  a 
constant  value.  In  this  limit,  the  radical  con¬ 
sumption  zone  retreats  far  upstream  into  the 
inert  so  that  the  fuel  consumption  region  is 
governed  solely  by  the  concentration  of  the 
fuel.  As  a  result,  a  coordinate  transformation 
can  be  made  which  expresses  the  structure 
problem  in  a  form  identical  to  that  found  by 
Lin4n  [9]  in  the  premixed  burning  regime  of 
his  diffusion  flame  analysis.  Thus  we  are  able 
to  adopt  his  correlation  curve  as  an  approxi¬ 
mate  expression  for  the  burning  rate  in  this 
limit.  We  find  that  the  role  of  the  propagation 
step  changes  from  retarding  to  enhancing  reac¬ 
tion  as  this  limit  is  approached. 

In  what  follows,  we  present  the  analysis  of 
the  above  structure  to  determine  the  charac¬ 
teristics  of  premixed  methane-air  flames.  In 
the  next  section  the  mathematical  formulation 
of  the  model  problem  will  be  discussed.  The 
asymptotic  analyses  of  all  the  sublayers  in  the 
flame  structure  are  then  performed  in  Sections 
3-5.  In  Section  6  we  analyze  the  resulting 
structure  equation  and  in  Section  7  we  evalu¬ 
ate  the  burning  rate  over  a  range  of  parameter 
values.  Finally,  in  Section  8  we  summarize  our 
results  and  add  further  discussion. 

2.  THE  MATHEMATICAL  MODEL 

The  reduced  four-step  mechanism  derived  in 
Ref.  3  consists  of  the  following  overall  reaction 


steps: 

CH4  +  2H  +  H,0  CO  -t-  4H,,  (I) 

CO -I- H,0  **  CO,  +  H,,  (II) 

2H  -t-  M  -*  H,  +  M,  (III) 

O,  +  3H,  2H,0  +  2H.  (IV) 


which  represent,  respectively,  fuel  consump¬ 
tion,  the  water-gas  shift,  radical  recombina¬ 
tion,  and  a  branching  reaction.  We  employ  the 


notation  of  Peters  and  Williams  [4]  to  express 
the  rates  of  these  four  steps  as 

w,  =  k„lCH  J[ri], 

(2.1a) 

A3li*2  J 

-[COjllHj/A;,), 

(2.1b) 

w„,  =  )tj[02][H][M]  =  k,„p[0^][H], 

(2.1c) 

(2.1d) 

where  p  is  the  pressure,  are  the  equilib¬ 
rium  constants,  and  kj  are  the  elementary  rate 
constants,  which  are  given  explicitly  in  Ref.  3. 
The  concentration  of  the  third  body  M  in 
reaction  III  has  been  written  in  terms  of  the 
Chaperon  efficiency  of  species  i  in  order  to 
convert  the  three-body  rate  constant  into  a 
two-body  form.  In  this  model,  reactions  I,  II, 
and  IV  htive  large  activation  energies  while 
reaction  III  has  zero  activation  energy  [3].  This 
implies  that  the  three-body  radical  recombina¬ 
tion  reaction  is  least  sensitive  to  temperature 
variations  and  therefore  takes  place  over  a 
broader,  but  still  asymptotically  thin,  region. 

We  start  out  by  first  assuming  that  the  wa¬ 
ter-gas  shift  (reaction  II)  is  in  partial  equilib¬ 
rium,  thereby  reducing  the  four-step  to  a 
three-step  mechanism.  This  assumption  can  be 
justified  when  the  corresponding  Damkdhier 
number  of  reaction  II  is  sufficiently  large,  a 
situation  that  can  be  achieved  at  sufficiently 
high  pressure.  We  consider  nonequilibrium 
effects  as  a  perturbation,  as  was  done  in  the 
studies  of  Peters  and  Williams  [4]  and  Seshadri 
and  Peters  [5].  That  is,  we  will  account  for  the 
existence  of  a  narrow  zone  immediately  down¬ 
stream  of  the  fuel  consumption  layer,  where 
the  water-gas  shift  is  not  in  equilibrium,  fol¬ 
lowed  by  a  layer  where  reaction  II  reaches 
equilibrium.  In  the  first  approximation,  we  set 
Eq.  2.1b  to  zero  to  obtain  an  expression  for  the 
concentration  of  CO. 

The  rate  constants  in  the  three  remaining 
global  reactions  govern  the  intermediate  reac- 
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tions 

CH4  +  H*-^CH3  +  H;, 

H  +  O,  ^  OH  +  O, 

H  +  O2  +  M  ^  HO,  +  M, 

which  represent  chain-propagating,  chain¬ 
branching  and  chain-terminating  steps  respec¬ 
tively.  Many  previous  studies  have  used  the 
rate-ratio  which  describes  the  competi¬ 

tion  of  oxygen  and  fuel  for  the  H  radical,  as 
the  small  perturbation  parameter  to  describe 
the  fuel  consumption  zone  thickness.  Those 
studies  have  found  that  ultimately  has  a 
chain-breaking  effect,  producing  the  more  sta¬ 
ble  methyl  radical,  CH  ,,  at  the  expense  of  the 
very  active  H  radical.  In  the  present  study,  we 
demonstrate  that  A:,,  can  play  an  intermediate 
role  to  that  of  the  chain-branching  and  chain¬ 
breaking  reactions  A:,  and  k,„,  respectively. 
Specifically,  we  find  that  when  a  particular 
lumped  parameter  containing  these  rates  is 
large,  the  chain-propagat  'g  reaction  further 
enhances  burning.  On  the  other  hand,  when 
this  lumped  parameter  is  small,  the  flame  bums 
weaker  and  reaction  A:,,  is  found  to  inhibit 
propagation. 

We  consider  steady,  planar  deflagrations  with 
constant  burning  rate  pi-  such  that  the  nondi- 
mensional  conservation  equations  for  the  re¬ 
maining  six  species,  as  well  as  the  temperature 
equation,  can  be  written  in  the  form 

Lr{Xf)^-w„  (2.2a) 

Z.^(A^^)  *  2o>y  2wm  2<sfjy,  (2.2b) 


In  addition,  the  partial  equilibrium  assumption 
for  reaction  II  yields  the  relation 

Xco  “ 

where 

^HjO^CO-^HjO*// 

Peters  and  Williams  [4]  argue  that  the  ratio 
^CO:  /A’hjO  varies  veiy  little  throughout  the 
flame  structure,  (insistent  with  their  study,  we 
treat  o  as  a  constant  in  order  to  facilitate  our 
analysis.  Here  the  differential  operators  are 
defined  as 


Li) 


d 


dP 


where  the  Lewis  number  of  each  species  Le, 
K/pDiCp  has  been  assumed  constant. 

ilie  nondimensional  quantities  in  our  ^- 
tem  are  given  in  terms  of  dimensional  vari¬ 
ables  (with  hats)  by 


and  the  reaction  terms  now  have  the  form 


Ly^^X^^,)  +  o^LQQiX■^^,) 

»  5a>,  ■¥  to,!,  -  3w,y, 


(2.2c) 


kWf 


I,..., IV. 


Lq^  Xq.)  =  tvyj-,  (2.2d) 

f"H;0^-^H.O^  “  ^L^Q(Xff^) 

■  -2<i>,  +  2iu,,  ,  (2.2e) 

Lco.^Xco.)  +  aLcoiX„.)  -  w,,  (2.20 

Ht)  +  Qii^Lf-QiX^^) 

“  <0/  +  0//)"/  +  Qiiit^ni  +  (2/1- w/i" 

.  (2.2g) 


In  the  above  equations  A  is  the  thermal  con¬ 
ductivity,  Cp  the  specific  heat,  and  pDt  the 
density-weighted  mass  diffusivities.  The  quanti¬ 
ties  Yi  and  If'  are,  respectively,  the  mass  frac¬ 
tion  and  the  molecular  weight  of  species  i  so 
that  each  Af,  is  related  to  the  mole  fraction 
normalized  with  respect  to  the  mole  fraction  of 
fuel  in  the  fresh  mixture.  The  subscript  u  de¬ 
notes  values  in  the  fresh  unbumed  mfarture 
and  f°  is  the  temperature  of  the  leading  edge 
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of  the  flame  zone  where  fuel  reacts  with  the  H 
radical.  The  parameter  Qi^  appearing  in  the 
energy  equation  (2.2g)  is  given  by 

^  ^F.uRk  .v 


CpW'rlfo  -  t)  ’ 


where  is  the  heat  release  per  mole  of  fuel 
consumed. 

The  above  equations  are  to  be  solved  subject 
to  the  boundary  conditions 

.Yf  =  Le;',  Jfo,  =  2Us;<^-',t=0,  1 

A’h  =  -^COj  “  -^HjO  *  “  ^CO  ““  1 

jf  _oe,  (2.5) 

■*^HjO  “  -^HjO.**  -^COj  ““  -^COj.*’ I 

A'o.  “  2Le3!(«5"'  -  1),  T  =  Tj,  I 
Xff  *=  Xf  =  •=  X^Q  =  0,  j 

z  -♦  *,  (2.6) 

where  d>  «  21^  „H'o./(yo,.„W'f)  is  the  equiva¬ 
lence  ratio,  such  that  d»  =‘  1  for  stoichiometric 
mixtures.  Our  analysis  can  only  be  considered 
valid  for  <^  <  1  since  the  reduced  mechanism 
and  the  corresponding  flame  structure  are  only 
appropriate  for  stoichiometric  and  fuel-lean 
mixtures.  An  analysis  of  fuel-rich  mixtures  will 
require  the  use  of  a  different  reaction  mecha¬ 
nism  and  flame  structure. 

It  is  convenient  to  work  with  the  reaction- 
free  coupling  functions 

L(t)  -t-  iQi  Q/i  Qiii  +  2Qiy)Lp{Xf) 


(Qiii  +  0/1') 


jaf-coC'^Hj^  “ 


L,(Xr)  +  K^(^//)  +  K„,(^„,) 

'*"5^HjO^^H;0^  *  (2.8) 

M-y,)  +  iL„(Ar„)-hiL„x^H:) 

“5^02('^o.^  ~  (2-9) 


which  replace  Eqs.  (2.2e,f,g).  These  equations 
can  be  integrated  across  the  entire  flow  field  to 
yield  the  following  leading  order  expressions 
for  the  burned  temperature  r^,  and  the  down¬ 
stream  product  concentrations  A’li^o.fr 

•^CO,.fc» 

Tj,  “  (2/  +  Qn  +  Qni  +  2Qiy,  (2.10) 

•^HjO.i  “  2LeH|o,  Xcoj.b  “  J-®coj* 

(2.11) 


Before  analyzing  the  ^tem  2.2,  we  first 
invoke  an  additional  assumption  in  order  to 
simplify  matters.  As  discussed  in  Ref.  4,  when 
the  H  radical  is  consumed  much  more  rapidly 
than  it  is  produced,  its  concentration  remains 
very  small  and  it  is  said  to  be  in  steady  state. 
Under  these  conditions,  the  above  system  of 
equations  is  simplified  and  the  three-step 
mechanism  is  further  reduced  to  a  two-step 
mechanism.  By  making  this  assumption  we  can 
set  the  right-hand  side  of  Eq.  2.2b  to  zero  to 
obtain 

lO/y  b>l  +  <Oiii,  (2.12) 

from  which  the  radical  concentration  is  found 
to  be 


Xff  “  jKjyXo^l^O^ 


LCh(-*HjO'^HjO 


^  ,  */// 
1^  ~  *,Uo,^02  ^ 


(2.13) 


In  Section  4  we  will  see  that  this  steady  state 
approximation  for  the  H  radical  holds  throu^- 
out  the  entire  flame  structure  with  the  mreep- 
tion  of  a  very  narrow  region  at  the  leading 
edge  where  ^e  fuel  attacks  the  radical  and 
totally  depletes  the  concentration  of  H. 

Equation  2.12  is  used  to  eliminate  w,y  so 
that  our  system  now  can  be  written 

Lf(Jff)--w,,  (2.14a) 

(Lhj  +  <*Lco)(-^Hj^  “  2(<tf;  —  '''//)* 

(2.14b) 

Lo/'^Oj^  “  ”**/  ~  (2.14c) 

^HjO^Xh^q)  —  aL^0(Xff^)  “  2<aiii,  (2.14d) 
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Lco,(^cO'^  +  “  <i>i,  (2.146) 

Ur)  -  Q„L»iX„^ 

“  (0;  ~  Qii  +  Qiy)<^i 

+  i2Q„  +  Q,„  +  Qn.)wi„.  (2.14f) 

The  remainder  of  this  study  is  concerned  with 
obtaining  solutions  to  Eqs.  2.14  in  order  to 
gain  a  description  of  the  methane-air  flame 
structure  sketched  in  Fig.  1;  we  begin  with  an 
analysis  of  the  chemically  inert  outer  regions. 

3.  THE  CHEMICALLY  INERT  REGIONS 

All  of  the  chemistry  is  confined  to  a  very 
narrow  region  located  near  j:  =  0,  on  either 
side  of  which  the  flow  field  is  chemically  inert. 
The  H  radicals  are  entirely  produced  and  con¬ 
sumed  within  the  flame  so  that  the  H-atom 
concentration  is  zero  on  either  side  of  the 
flame  zone.  For  the  remaining  species  a  bal¬ 
ance  is  maintained  between  convection  and 
diffusion.  Solutions  are  readily  found  to  be 

Le^ '(I  -  exp(Lefjr  )1  ^<0,  ) 

0  jr  >  0, 


layer,  respectively,  and  they  will  be  defined  in 
the  next  two  sections. 

We  note  that  the  fuel  is  entirely  consumed 
by  the  flame  so  that  no  leakage  to  the  down¬ 
stream  region  occurs.  Furthermore,  we  will 
impose  continuity  of  Xf  across  the  fuel  con¬ 
sumption  zone  up  to  0(£),  which  is  the  needed 
order  of  accuracy  for  the  present  problem. 
When  continuity  is  applied,  the  concentration 
of  fuel  outside  the  flame  structure  is  given 
entirely  by  Eq.  3.1  to  all  orders  in  €  and  S. 

4.  THE  FUEL  AND  RADICAL 
CONSUMPTION  LAYERS 

In  seeking  solutions  to  the  system  2.14  we 
assume  that  the  nonzero  activation  energies 
are  sufficiently  large  so  as  to  confine  reaction 
of  fuel  to  a  narrow  region  in  the  flow  field, 
thus  permitting  the  use  of  activation  energy 
asymptotics.  Of  the  reaction  rates  remaining  in 
our  system,  ki„  has  zero  activation  energy 
while  and  k,  have  the  form 

^11  "  ^ll.eff 

k, -/l,exp(-£,/R'>f),  (4.1) 


2Le3l<^' '  -  6,  exp(Leo,z)  x  <  0, 

2U3.'(<f."' -  1)  x>0. 


(>2  expd-en.Ar)  x  <  0, 

0  z  >  0, 

|f»3exp(LeH,o-f)  x<0, 

Jf  >  0, 


(3.2) 


(>4  exp(Leco,Jt)  x  <  0, 
^CO:.b  X  >  0, 


X  <  0, 

I  Tj  Soj  +  "•  ,  X  >  0, 


(3.3) 

(3.4) 


where  the  unknown  constants  a,,  6,  will  be 
determined  by'  matching  to  solutions  in  the 
inner  flame  region.  In  general,  all  constants 
are  to  be  expanded  in  terms  of  the  two  param¬ 
eters,  5  and  c,  which  are  assumed  small  rela¬ 
tive  to  the  thermal  thickness  of  the  flame. 
These  parameters  represent  the  thickness  of 
the  fuel  consumption  zone  and  the  oxidation 


where  the  E^s  are  the  activation  energies  and 
R°  is  the  gas  constant.  The  preexponentials 
are  generally  quite  large  and  so  it  is  appropri¬ 
ate  to  rescale  them  as 


A,^A,exp{Ei/R°f°),  (4.2) 

where  7°  was  defined  earlier  as  the  character¬ 
istic  temperature  at  which  fuel  reacts  with  the 
radical.  It  follows  from  Eq.  2.14a  that  fuel  is 
consumed  by  the  reaction  with  the  rate  k,, 
which  suggests  that  an  appropriate  small  pa¬ 
rameter  to  define  the  thickness  of  the  ffiel 
consumption  layer  is  given  by 


|jO(fO)2 

E{f°  -  f„)  ■ 


(4J) 


Thus,  fuel  consumption  will  be  restricted  to  a 
narrow  region  of  width  0(5)  where  the  tem¬ 
perature  is  close  to  f®,  that  is,  where  t  -  1  - 
0(5).  We  note  that,  although  the  temperature 
will  continue  to  increase  downstream  of  this 
layer,  reaction  is  terminated  by  the  depletion 
of  fuel. 
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4.1.  The  Fuel  Consumption  Layer 

In  terms  of  the  outer  variable,  the  flame  re* 
sides  at  the  location  jr  0.  To  investigate  the 
fuel  consumption  region  we  introduce  the 
stretched  coordinate  x  •=  S(  and  seek  solu¬ 
tions  of  the  form 

T  »=  1  +  ,  A"/:  •=  Syf  ••• , 

A,  *  A°  +  Sy,  +  • ,  (i  #  F).  (4.4) 

The  H-radical  concentration  is  readily  deter¬ 
mined  from  Eq.  2.13  as 


fuel  consumption  zone  (2.14a,  b)  and  (2.7)  can 
be  expressed  as 


(4.6) 


+  RfVf  +  9h>’h 


+  (9h,  +  «9co)>’H:]  *=  0. 
d^{\  +  a)yH,  ^d^yF 
dC^' 


(4.7) 

(4.8) 


where  we  have  defined  the  parameters 


, - - -  (LeH  A^  ) 


^iI-®o.Ao,  ^1 


.  (4.5) 


Gradients  are  large  in  this  region  so  that  dif¬ 
fusion  dominates  over  convection  and  thus  we 
expect  the  diffusion  and  reaction  terms  to  bal¬ 
ance  in  our  equations.  In  order  that  such  a 
balance  is  maintained  in  the  equation  for  the 
temperature  perturbation,  we  require  that  A;,, 
*  0(6”^)  so  that  A,,  “  We  also  want 

to  retain  the  variation  of  H  radical  with  fuel 
concentration  and  by  balancing  the  first  two 
terms  on  the  right-hand  side  of  Eq.  4.5  it  is 
necessary  to  introduce  the  additional  scaling 
A,  *  5" 'A,.  Furthermore,  the  ratio  A:,//c,]  has 
a  weak  temperature  dependence  [4].  Thus  we 
simplify  matters  by  assuming  the  activation 
energies  of  these  two  rates  are  equal,  that  is, 
Ej)  »£,*£.  The  above  scalings  imply  that 
the  last  term  in  Eq.  4.5  is  smaller  in  magnitude 
than  the  remaining  terms  throughout  this  re¬ 
gion  (i.e.,  k,„/k^  <«  1),  and  that  the  ratio 
k\/kji  ^  6A^/A^J  is  a  small  parameter.  In¬ 
deed,  as  discussed  earlier,  many  previous  stud¬ 
ies  have  used  the  rate-ratio  A:, /A:,,,  which  de¬ 
scribes  the  competition  of  oxygen  and  fuel  for 
the  H  radical,  as  the  small  perturbation  param¬ 
eter  to  describe  the  fuel  consumption  zone 
thickness. 

To  leading  order  in  5,  reaction  w,,/  is  negli¬ 
gible  and  the  relevant  equations  governing  the 


A  =  Ai]Le^/(AiLeojAoj), 

UfAy>_„p^  r 


p^v^cWf 


lK,yX^:Uo, 


^  1-  vo 

^HjO-f'HiO 

"  G/  +  Qi,  +  Qm  +  20,,., 

9h  “  ^^Qin  G/kV^* 

9h,  “  ^G///  Qiv^/^^ 

9co  “  G;/  +  (G///  G/k)/2. 


(4.9) 


(4.10) 


These  equations  are  to  be  solved  subject  to 
appropriate  boundary  conditions  which  are  ob¬ 
tained  through  matching  with  the  upstream 
inert  region  and  the  downstream  oxidation 
layer.  We  have  avoided  writing  equations  for 
the  Oj,  HjO,  and  COj  perturbations  since  the 
system  4.6-4.8  is  independent  of  these  quanti¬ 
ties. 

As  we  mentioned  earlier,  the  steady-state 
approximation  implies  that  the  H-radical  con¬ 
centration  is  small  so  that  can  be  neglected 
in  Eq.  4.7.  However,  it  is  easily  seen  that  Eq. 
4.6  breaks  down  as  -»  1/A,  suggesting  that 
our  steady-state  assumption  for  H  becomes 
invalid,  llius  we  expect  a  boundary  layer  to 
develop  near  y^  “  1/A,  where  the  fuel  de¬ 
pletes  the  radical  concentration.  Peters  and 
Williams  [4]  have  analyzed  this  boundary  layer 
and  have  demonstrated  the  smooth  transition 
to  a  zero  H-atom  concentration.  For  complete¬ 
ness,  we  briefly  summarize  that  analysis. 
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A2.  The  Radical  Consumption  Layer 

We  first  need  to  determine  the  local  behavior 
of  Xf.  Since  the  region  preceding  the  radical 
consumption  layer  is  chemically  inert  (due  to 
the  absence  of  the  H-radical),  the  solution  3.1 
for  Xf  is  valid  everywhere  upstream  of  the 
boundary  layer  at  ^  -\ /A.  TIjus  by  expand¬ 

ing  Eq.  3.1  in  terms  of  the  inner  variable 
jr  “  we  find 

(4.11) 

from  which  it  follows 


for  -1/A.  (4.12) 

In  the  remainder  of  the  structure,  i.e..  (  > 
-1/A,  y'f  is  determined  by  solving  Eq.  4.6, 
and  continuity  of  »  and  d\f  /d^  at  f  =  -1/A 
provides  the  boundary  conditions 

ziC~-l/A.  (4.13) 


The  equation  governing  the  radical  concen¬ 
tration  in  the  flame  structure  (where  a;,,,  is 
negligible)  (Eq.  2.2b)  can  be  written  as 


d^n 

d(- 


2d~(.u},  -  a>,,.) 


s- 


1  ^LCq.Lch  ATf  „  Pjl 


p-i-CfWf 


k  y®  y® 


X 


-1  -I- 


k^^LCfXfr 

^ll-*0;'^0; 


+  (^H^H:0-^H:0-^h) 

LeQ^.ArO,/(/,. 

(4.14) 


In  the  much  thinner  radical  consumption  layer 
temperature  variations  are  negligible,  so  that 
6  »  0,  and  we  introduce  the  stretched  variable 
o  as 

-1/A  +  ^Jia/A,  (4.15) 

where  the  small  parameter  p  will  be  appropri¬ 
ately  chosen  shortly.  It  follows  frbm  Eq.  4.13 


that  the  fuel  concentration  in  this  layer  has  the 
form  yf  -  l/>4  -  fta/A.  We  insert  this  into 
the  steady-state  expression  for  (4.5)  to 
obtain  the  following  matching  condition  for 


A®  -- 


^ KivLcq^Xq^  (LeHj-^Hj) 
^hLChjO-^HjO 


O’  -»  «.  (4.16) 

At  the  opposite  end  of  this  region  the  radical  is 
attacked  and  consumed  by  the  fuel,  yielding 
the  additional  boundary  condition 


AS-0,  O’-*-*. 


(4.17) 


The  matching  condition  (4.16)  suggests  the 
magnitude  of  to  be  OfyfiL)  in  this  region 
and  it  is  convenient  to  introduce  the  new  vari¬ 
able  W  as 


AS 


•^KiylxQ^Xo,  (LChjAh,) 
LChLCh.O-^HjO 


(4.18) 


so  that  our  system  governing  the  radical  con¬ 
sumption  layer  takes  the  simple  form 


d^W 

0 


W[W-  -  O’], 

as  o  -* 


(4.19) 

IF  =  r/a  as  O’  -♦  «. 

(4.20) 


Here  p  is  chosen  as 

I  A^  ^A/,/LeojAoj  (Lch^Ahj)  | 

^  |26D  j 

«  1,  (4.21) 

where  we  have  made  use  of  definition  4.10. 
The  system  4.19-4.20  can  be  solved  numeri¬ 
cally  [4]  and  the  resulting  profiles  for  fV  illus¬ 
trate  the  smooth  transition  to  zero  of  the  H- 
atom  concentration  in  the  radical  consumption 
zone.  . 

We  remark  that  the  analysis  of  this  region  is 
based  on  the  assumption  that  p  is  a  small 
parameter  and  it  follows  from  Eq.  4.21  that 

Leo.Ao,  ^LeH.AH,) 

LcHLeH.o-^H.o 


S.  (4.22) 
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By  comparing  this  expression  to  Eq.  4.5  we 
conclude  that  the  radical  concentration  must 
be  smaller  in  magnitude  than  the  width  of  the 
fuel  consumption  zone,  which  is  in  agreement 
with  our  earlier  comments  regarding  the  valid¬ 
ity  of  the  steady-state  approximation. 

As  discussed  in  Ref.  4,  the  existence  of  this 
radical  consumption  layer  indicates  that  it  is 
the  consumption  of  H  atoms  by  the  fuel  that 
terminates  the  reaction  in  the  upstream  pre¬ 
heat  zone.  This  differs  from  standard  analyses 
with  one-step  kinetics  where  reaction  is  frozen 
in  the  preheat  zone  by  a  drop  in  temperature. 

Based  on  the  above  boundary  layer  analysis 
we  can  neglect  the  perturbation  in  Eq.  4.7. 
Furthermore,  we  can  use  Eq.  4.8  to  replace 
and  thus  reduce  the  structure  of  the  fuel 
consumption  layer  to  the  two  equations 

, - 

-jp-  “  Dy-fyjl  -  Ay'f  exp(e),  (4.23) 

£_ 

d(‘ 

(4.24) 


+  a<lcof 

0  +  9r>r  -  2 - — - » 


=  0. 


Upstream  boundary  conditions  for  this  system 
are  obtained  by  matching  to  the  solutions  in 
the  preheat  zone,  where  in  particular  we  have 
obtained  Eq.  4.13  for  ,  and  the  expansion  of 
the  outer  solution  3.4  in  terms  of  x  -  5^  pro¬ 
vides  the  necessary'  conditions  for  6.  The  fuel 
is  entirely  consumed  downstream,  which  im¬ 
plies  that 

»  -  0  as  -♦  «,  (4.25) 

and  boundary  conditions  for  6  are  obtained  by 
matching  to  the  solution  in  the  oxidation  layer 
immediately  adjacent  to  the  fuel  consumption 
zone.  An  analysis  of  this  oxidation  layer,  which 
we  now  present,  also  provides  us  with  the 
needed  expression  for  appearing  in  D,  as 
shown  in  Eq.  4.10. 


which  will  be  chosen  shortly  such  that  See 
‘C  1.  For  this  reason  the  present  flame  struc¬ 
ture  resembles  that  of  the  Zerdovich-Linan 
mechanism  [10,  11],  which  consists  of  a  very 
thin  branching  reaction  zone,  characterized 
a  large  activation  energy,  embedded  within  a 
somewhat  broader  recombination  zone*  with 
zero  activation  energy. 

Our  governing  equations  simplify  in  this  re¬ 
gion  since  -  0  and  ta,  ~  0.  The  molar  con¬ 
centrations  for  H  and  CO  are  obtained  directly 
from  Eqs.  2.13  and  2.3.  In  order  to  determine 
the  remaining  species  concentrations  we  intro¬ 
duce  the  stretched  coordinate,  17,  as 


X  =  CT)/2q, 

(5.1) 

and  seek  solutions  of  the  form 

T=  Tfc  -  6/  +  •", 

q{\  -f  a)A’H,  =  cZ  -1-  •••, 

•^Oj  ~  ^Oi.b 

(5.2) 

■^H.O  “  ■^HJO.^  ”  > 

^COi  “  -^COj.fc  ~  ®2cO:  ■■■  • 

(5.3) 

At  leading  order,  the  equations  governing  this 
region  can  now  be  written  in  the  form 


d^Z  _  € 
dr)^  2q 


d‘ 


2  -  ^[Z  -  2qZo.]  =  0, 


d^ 

-^[Z-t]  =  0, 

dri^ 

d-  ,  ^  , 

^^|Z  -  q(l  +  a)ZH,o]  * 


il 

drf^ 


(1  +  a) 

Z  -  q  -  Zco. 


0, 


(5.4) 


where  we  have  employed  the  notation 


5.  THE  Hj  AND  CO  OXIDATION  LAYER 

Immediately  downstream  of  the  fuel  consump¬ 
tion  region  lies  a  somewhat  broader  layer  where 
H,  and  CO  oxidize  to  form  H,0  and  COj. 
The  width  of  this  layer  is  characterized  by  e 


9  “  9hi  “  9co  “  ^Qm  Qiv^/^-  (5.5) 

Peters  and  Williams  [4]  have  computed  the 
values  of  found  that  Q,,  is  sufficiently 

small,  thus  enabling  us  to  write  the  second 
equality  in  Eq.  5.5.  By  integrating  these  equa- 
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tions  and  matching  to  the  outer  inert  solutions 
3. 1-3.4  at  rt  -*»  we  obtain  the  relations 

S 

t  “  Z  —  Zq^ 

Z 

“  q(l  +  a)' 


where  a  =  4q(l  -  <^)/(«LeQ^)  for  near- 
stoichiometric  flames.  These  are  all  ejqiressed 
in  terms  of  the  H;  concentration,  2,  which  in 
turn  is  found  by  solving  Eq.  5.4a.  The  nature  of 
this  structure  equation  changes  depending  on 
whether  the  mixture  is  near-stoichiometric  or 
if  there  is  an  abundance  of  oxygen.  The  analy¬ 
sis  of  this  region  for  stoichiometric  flames  has 
been  performed  by  Peters  and  Williams  [4], 
and  a  similar  analysis  for  near-stoichiometric 
flames  can  be  found  in  the  analysis  of  Seshadri 
and  Gottgens  [6]. 

5.1.  The  Oxidation  Layer  for  Near- 
Stoichiometric  Mixtures 


2  +  n 

‘~W’ 

qZ 

“  q(Wa)’ 

(5.6) 


where  e  has  been  chosen  to  be 


P^ni 


AyliqW  4- 


-J/4 


(5.10) 


Alternatively,  we  can  use  Eq.  4.10  to  express  e 
in  terms  of  0(1)  quantities,  which  yields 


j  4\/2q^ 

I  DLCqj  pkn, 


Since  e  must  be  a  small  parameter,  it  follows 
that  pk,„  must  be  sufficiently  large.  We  re¬ 
mark,  however,  that  it  must  not  be  so  large  as 
to  violate  our  previous  assumption  that  k,„  c 
k,  -  0(6-'). 

Boundary  conditions  for  Eq.  5.9  arc  ob¬ 
tained  by  matching  to  solutions  in  the  outer 
inert  region  at  t)  and  to  the  solutions  in 
the  fuel  consumption  layer  at  17 »  0.  This 
matching  procedure  yields 


When  1  -  d>  “  0(«),  the  oxygen  is  nearly 
entirely  consumed  by  the  flame,  so  that  the 
oxygen  concentration  is  at  most  0(e)  in  the 
oxidation  layer.  We  use  Eq.  2.1c  to  write  the 
equation  for  2  in  the  form 


dZ 

—  *  - 1  at  17  ■=  0, 

dv 

(5.12) 

dZ 

—  “2  =  0  asi7-»*, 
di) 

(5.13) 

d^Z  e  Lco.LCh  ATf  „  P5  ^ 
dv-  "  2q 

(5.7) 


and  the  H-radical  concentration  obtained  from 
Eq.  2.13  is  given  by 


^X/i/LCq,  Le  i(‘ 

^H^HjO-^HjO-^H:0,fc 


}/Xn'l^o,LeK^  e^Z^^HZ  + 

^hLchjO-^HjO.s  'Ilq^d  + 

(5.8) 


This  is  now  inserted  into  Eq.  5.7  to  obtain  the 
equation 


d^Z 


Z^''(Z  -t-  a)^^\ 


(5.9) 


where  the  first  of  these  is  obtained  by  a  single 
integration  of  Eq.  4.8  across  the  entire  fuel 
consumption  region.  For  the  stoichiometric 
flame,  a  >  0,  the  solution  of  Eq.  5.9  subject  to 
these  boundary  conditions  is  readily  found  to 
be 


n 

^  “  1,  -H  2'/“  ’ 


(5.14) 


from  which  the  concentration  in  the  fuel 
consumption  zone  is  determined  from  2(0)  -> 
2'''*.  For  nonzero  values  of  a,  a  single  integra¬ 
tion  of  the  above  system  determines  2(0)  as 
2(0)  -  fl(j  -  l)/2,  where  s  is  the  root  of  the 
equation 

s}/s^  -  1  +  |ln[s  -I-  -  1  ]  -  32/a". 

Solutions  for  the  remaining  quantities  follow 
immediately  from  Eqs.  5.6,  and  their  evalua- 
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lion  at  17  »=  0  provides  the  necessary  matching 
conditions  for  variables  in  the  fuel  consump¬ 
tion  layer. 

The  above  analysis  assumes  partial  equilib¬ 
rium  of  the  water-gas  shift.  However,  as  dis¬ 
cussed  in  Ref.  4,  assuming  this  reaction  to  be 
infinitely  fast  overestimates  the  amount  of  Hj 
produced.  A  correction  due  to  nonequilibrium 
effects  can  be  obtained  by  considering  an  addi¬ 
tional  layer  of  Ofi"),  i*  •«  e  «  1,  on  the  down¬ 
stream  side  of  the  fuel  consumption  zone  where 
reaction  II  is  not  in  equilibrium.  Peters  and 
Williams  [4]  have  performed  that  analysis  so 
the  details  will  not  be  presented  here.  To  ex¬ 
amine  this  boundary  layer  one  first  perturbs 
the  mass  fractions  of  H,  and  CO,  which  are 
the  relevant  quantities  in  this  zone,  about  their 
value  at  Tj  =  0.  The  resulting  equations  for  the 
perturbed  quantities  are  linear  and  can  easily 
be  integrated  subject  to  appropriate  matching 
conditions  [4].  This  procedure  is  straightfor¬ 
ward  and,  when  the  solution  is  combined  with 
that  of  the  partial  equilibrium  layer,  the  con¬ 
centration  of  H2  at  the  fuel  consumption  zone 
is  found  to  be 


eZ(0) 
q(l  +  a) 


t'. 


(5.15) 


the  flame  unreacted.  The  resulting  flame  struc¬ 
ture  has  not  been  considered  in  previous 
asymptotic  studies  that  have  employed  rate- 
ratio  techniques.  In  this  case,  the  structure 
equation  5.9  is  modified  since,  to  a  first  approx¬ 
imation,  the  molar  concentration  .Yq- 
the  value  j,  =  0(1),  shown  in  Eq.  3.2. 

All  the  relations  in  Eqs.  5.6  remain  un¬ 
changed,  but  now  the  leading  order  expression 
for  the  radical  concentration  is  given  by 


iJKiv'^oMh^  ,/2 

ije.  y  '■■^o,.*' 


(5.16) 


and  when  this  is  inserted  into  Eq.  5.7  we  ob¬ 
tain 


d-Z  €  Le^fAT^.^pg 

dr,^  2q  ph-^CpWf  Lch^o^HjO.* 


cZ 

q(\  -I-  a)  / 


(5.17) 


where 


1  —  0 
(1  +  o)’  '- 


(1  +  a)‘  1 9  / 

^  j  LcQ.LeH.M///^?  \ 

\  8UcoA,oVZ(0)[Z(0)  +  a]  j 

The  expression  5.15  provides  the  matching 
condition  for  the  variables  in  the  fuel  con¬ 
sumption  layer  and  will  be  applied  in  Section  6, 
but  we  first  reexamine  the  oxidation  layer  for 
the  case  when  the  mixture  contains  an  abun¬ 
dance  of  oxygen. 


The  boundary  conditions  5.12-5.13  are  un¬ 
changed  so  that  our  system  takes  the  simple 
form 


d-Z 

drj- 


at  Tj 


0, 


-;-«=Z  =  0  aST}->=c, 

dv 

where  e  is  now  given  by 

I  .11  PO^^O;^H;^  '' 

I  LCh.-O-^HjO.* 


(5.18) 

(5.19) 

(5.20) 


5.2.  The  Oxidation  Layer  for 
Fuel-Lean  Mixtures 

When  d>  <  1  there  is  a  surplus  of  oxygen  in  the 
mixture  so  that  a  portion  of  it  passes  through 


pk 


III 


2q’/-(l  -t-  a) 


3/: 


3/2 


-2/5 


«  1, 


(5.21) 
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dr  alternatively 

' - 7^-  (5.22) 

2JLeo,Ao,.j,  pk/ii 

As  in  the  case  of  the  stoichiometric  mixture, 
the  requirement  that  c  be  small  implies  that 
the  three-body  reaction  rate  pk,,i  is  large, 
since  all  other  quantities  in  Eq.  5.22  are  known 
to  be  0(1). 

The  system  5.18-5.20  possesses  the  solution 


(5.23) 


from  which  we  obtain  the  needed  information 
at  t;  =  0,  Z(0)  =  (v/5/2)'''-. 

As  in  the  preceding  subsection,  we  can  con¬ 
sider  nonequilibrium  effects  of  reaction  II  as  a 
perturbation  to  the  above  leading  order  result. 
Again  Eq.  5.15  is  obtained  for  the  H;  concen¬ 
tration  at  the  fuel  consumption  zone  where 
now 


given  entirely  by  Eqs.  3. 1-3.4.  Downstream  of 
the  fuel  consumption  layer,  the  solutions  in  the 
inert  region  and  the  oxidation  layer  can  be 
combined  to  construct  composite  solutions  that 
are  valid  in  the  entire  flow  field.  We  now 
discuss  these  composite  solutions,  as  well  as 
the  resulting  matching  conditions,  for  both  sto¬ 
ichiometric  and  fuel-lean  mixtures. 


6.1.  Structure  Equatlcn  for  Near- 
*  Stoichiometric  Mixtures 

For  near-stoichiometric  mbctures,  the  solutions 
downstream  of  the  fuel  consumption  layer  are 
given  by  Eqs.  3.1-3.4,  5.6,  and  5.14.  Appropri¬ 
ate  composite  solutions  for  t*,  Xf  and 
are  given  by 


x;  s  0, 


X 


H. 


Tj,  +  So,  -  cZ  -I-  ••• , 

'  '^oTT)  ^  ■■■■ 


(6.1) 

(6.2) 

(6.3) 


1  —  Q 

(1  -I-  af  '-y/Di; 

1  -  a  € 

(1  +  o)^'-  q 

^\4(1  -t-  o)Leeok,oVZ(0)  j 


To  leading  order  in  S  we  impose  continuity 
at  X  0,  and  the  nondimensionalization  fur¬ 
ther  requires  that  t(0*)  =  t(0“)  =  1,  which 
yields 


a,  =  1, 


6,  =  e 


Z(0) 

q(l  +  a)’ 


Tf  =  1  +  eZ(0)  -I-  o(c). 


(6.4) 


We  now  have  the  necessary  matching  condi¬ 
tions  to  complete  the  analysis  of  the  fuel  con¬ 
sumption  region  that  was  begun  in  Section  4. 

6.  THE  FINAL  STRUCTURE  EQUATION 

In  Sections  3  and  5  we  obtained  solutions  for 
each  species  as  well  as  the  temperature  distri¬ 
bution  in  both  the  inert  regions  and  the  H  i-CO 
oxidation  layer.  These  solutions  provide  the 
matching  conditions  needed  to  solve  the  struc¬ 
ture  equations  4.23  and  4.24.  For  our  purpose, 
we  only  need  to  consider  the  expressions  for  t, 
Xf  and  A'h..  Upstream  of  the  fuel  consump¬ 
tion  layer,  solutions  for  these  quantities  are 


The  last  of  these  relations,  together  with  Eq. 
2.10,  provides  the  two  equations  necessary  to 
determine  the  dimensional  temperatures  f|, 
and  f 

The  matching  conditions  needed  for  Eqs. 
4.23-4.24  are  now  obtained  by  expanding  the 
outer  solutions  in  terms  of  the  inner  variable 
X  -  6{.  At  this  stage  of  the  analysis,  we  have 
not  explicitly  related  the  two  independent  small 
parameters,  e  and  5,  although  we  have  chosen 
them  such  that  6  •«  e.  The  matching  proce¬ 
dure  requires  that  we  retain  solutions  up  to 
0(5^),  and  thus  we  will  now  assume  that  the 
solutions  5.6  are  valid  to  all  orders  in  e.  We 
neglect  all  terms  of  o(6^)  so  that  the  matching 
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conditions  are 

X}  -  6V'Lef/2+  •••, 

+  52^'LeJ,y2 +  •••), 


T~“»  1  +  SCf  +  flj)  +  S~  +  "■ » 


-  0  (6.8) 

+  6'^29'2Z,,(0)/c+ -).  (6.9) 

T*-  1  +  Sla^  +  2qn  -2S-iYZ^r,(0)/e] 

+  -.  (6.10) 

Integration  of  Eqs.  4.23  and  4.24  subject  to  the 
above  ntatching  conditions  now  yields  the  lead¬ 
ing  order  result 

e~  -yF(qr~2q)  +  2qi,  (6.11) 

where  02  *  flj  =  0.  We  insert  Eq.  6.11  into  Eq. 
4.23  to  obtain  the  following  boundary  value 
problem  for  the  fuel  concentration  in  the  fuel 
consumption  layer 

d^y'F  I - 

-jp-  »  Dyry/1  -  Ay, 

xexpl  -»(9f  -  2q)  -t-  2q{  ], 

(6.12) 


y^  -  -i-  at  -1/A. 

»  “  0  as  f  -*  X.  (6.13) 

The  burning  rate  pv  appears  as  an  eigenvalue 
of  this  system.  In  particular,  for  stoichiometric 
mixtures  we  note  from  Eq.  5.10  the  relation 
e  ~  i/pT,  and  it  follows  from  Eqs.  4.9  and  5.11 
that  the  rate-ratio  parameters  A  and  D  de¬ 
pend  on  the  burning  rate  as 

A  1/v^,  D  '  l/pf. 


In  general.  Eqs.  6.12  and  6.13  (which  includes 
the  evaluation  of  the  burning  rate  in  terms  of 
all  remaining  parameters)  must  be  solved  nu¬ 
merically.  However,  we  will  firrt  discuss  several 
limiting  cases  for  which  explicit  expressions  for 


the  burning  rate  can  be  obtained.  Computed 
values  of  .  the  burning  rate  as  a  function  of 
various  parameters  will  be  given  in  Section  7. 

We  first  consider  the  limit  A  -c  1.  This  can 
be  achieved  for  example  when  the  chain¬ 
branching  reaction  becomes  comparable  to  the 
chain-propagating  step  such  that  A', 

We  should  mention  that  such  small  values  of 
A  are  unlikely  for  the  methane-air  flames  of 
interest  here,  and  the  steady-state  assumption 
for  the  H  radical  becomes  questionable.  How¬ 
ever,  it  is  possible  to  achieve  smaller  values  of 
A  by  using  a  different  inert  or  by  considering  a 
different  hydrocarbon,  for  which  the  basic 
flame  structure  remains  the  same  (cf.  the  anal¬ 
ysis  of  propane-air  flames  in  Ref.  12).  In  any 
case,  it  is  instructive  to  examine  this  limit  in 
order  to  determine  the  explicit  dependence  of 
the  burning  rate  on  various  parameters.  In  this 
limit,  the  radical  consumption  layer  moves  far 
upstream  into  the  preheat  zone;  the  radical 
concentration  assumes  the  value  unity  (i.e., 
y/i  -  Ay f  1);  and  the  fuel  consumption  layer 
is  limited  by  a  single  reactant  namely  the  fuel. 
In  fact,  we  can  make  the  coordinate  transfor¬ 
mations 


/f  “  yA<iF  -  29), 

9f  -  2q 

1  .  [  -  29)^ 

+  —In  - — - 

2q  2D 


(6.14) 


and  introduce  the  parameter  m  -  2q/(2q  — 
qf)  to  write  our  system  in  the  form 


d^Y 

2-jp  •=  Texpl-y-  ], 


(6.15) 


1  29* 


as  ^  -<*, 


y  -  0  as  f  -» X, 


(6.16) 


which  is  identical  to  the  structure  problem 
obtained  by  Lihan  [9]  in  the  premixed  burning 
regime  of  his  diffusion  flame  analysis.  The 
latter  boundary  condition  implies  there  can  be 
no  leakage  of  fuel  through  the  flame.  Ijndn 
has  demonstrated  that  solutions  to  the  above 
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^stem  only  exist  when  -*  <  m  <  0.  In 
methane-air  flames,  the  value  of  m  can  be 
taken  as  m  ~  -  2  [4],  and  thus  we  are  in  the 
range  of  physically  meaningful  solutions.  Linan 
[9]  has  also  provided  a  correlation  function  for 
the  limiting  value  lim^__,(y  +  ^)  which  we 
use  to  obtain  the  following  expression  for  D 

2q^ 

D  -  -^|0.6307m2  -  1.344m  +  1].  (6.17) 

m* 


where 

-TH/7-m  +  2m/), 

/  s  J  wCl  —  w)!!  +  2w  +  3^^  +  dw 

-  0.291...  , 

From  the  definitions  of  D  and  A,  we  find  that 
the  rate-ratio  parameter  D  has  the  e}q>licit 
form 


Now  the  definition  of  D,  Eq.  4.10,  can  be  used 
in  this  equation  to  obtain  an  expression  for  the 
burning  rate,  pv,  for  small  values  of  A.  We 
note  that  since  m  <  0,  the  solutions  of  the 
above  ^stem  do  not  exhibit  turning  point  be¬ 
havior,  and  thus  no  extinction  is  predicted. 

In  the  opposite  limit  that  A  radical 
production  is  extremely  weak,  and  thus  we 
expect  the  flame  to  slow  down  and  possibly 
“extinguish.”  To  analyze  this  limit  it  is  conve¬ 
nient  to  remove  the  parameters  from  the 
boundary  conditions  by  introducing  the  scaled 
variables 

Xf^y/A.  (6.18) 

Our  system  can  then  be  written 
d*  * 

jp  -  Dy y/l  -yexp[-y(y  -(■  m/)].  (6.19) 

>•  -  -  /  at  /  =  - 1,  y  «  0  as  /  -♦  X, 

(6.20) 


where  we  have  defined  the  parameters 

9f  -  2q  2q 

ym  - - — ,  mm- - , 

A  2q  -qf 

D  m  D/A-.  (6.21) 

We  can  now  exploit  the  limit  of  small  y. 
When  y  is  sefto  zero  in  the  above  system  we 
recover  the  structure  problem  of  Peters  and 
Williams  [4].  The  resulting  equation  can  be 
integrated  once  subject  to  the  boundary  condi¬ 
tions  to  yield  D  *  15/8.  A  correction  for  small 
values  of  y  is  readily  found  to  be 

D'15/8+yD,.  .  (6.22) 


D 


A^Uq,  Z(0f''^(Z(0)  + 
P^III 


1  - 


vq{\  +  a) 
€Z(0) 


3/2 


(6.23) 


If  we  consider  stoichiometric  flames  (a  »  0) 
and  invoke  the  partial  equilibrium  assumption 
(I'-O),  then  the  leading  order  result  D  •> 
15/8  does  not  provide  an  expression  for  the 
burning  rate.  Therefore  Peters  and  Williams 
[4]  used  this  leading  order  result  to  determine 
f  and  then  equated  the  first  two  terms  in  the 
expansion  for  r^,,  shown  in  Eq.  6.4,  to  evaluate 
the  burning  velocity.  In  our  model,  this  value 
of  D  is  approached  in  the  limit  that  y  -*  0. 
Since  y  is  proportional  to  the  square  root  of 
the  burning  rate  (see  Eq.  6.21),  D  «  15/8 
represents  the  critical  point  at  which  pv  de¬ 
creases  to  zero.  In  this  limit  e  -*  0,  the  oxida¬ 
tion  layer  collapses  down  to  the  fuel  consump¬ 
tion  layer,  and  the  burned  temperature,  r^, 
decreases  to  unity. 

In  Fig.  2  we  have  plotted  y,  which  repre¬ 
sents  the  burning  rate  for  the  stoichiometric 


Fig.  2.  Plot  of  y,  representing  the  burning  rate  for  noi- 
chiometric  flames,  versus  the  rate-ratio,  D.  Curve  is  gener¬ 
ated  b)'  soK'ing  equations  6.19  and  6.20. 
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flame,  as  a  function  of  D  for  the  fixed  value 
m  •=  -2.  It  is  clear  from  Eq.  6.23  that  the 
parameter  D  represents  the  ratio  of  the 
chain-branching  reaction,  A:,,  to  the  chain- 
propagating  and  chain-breaking  reactions,  k,, 
and  k,,,.  As  previously  discussed,  when  D  is 
sufficiently  large  the  flame  has  a  structure  sim¬ 
ilar  to  the  premixed  flame  regime  of  Lilian’s  [9] 
diffusion  flame.  The  limiting  value  of  the  burn¬ 
ing  rate  is  obtained  from  Eq.  (6.17),  which  can 
Ik  expressed  in  terms  of  the  parameters  y  and 
D  as 

,  ,  -1  2D 

■y*  =  2D[0.6307m*  -  1.344m  -t-  1]  =  - . 

M 

For  the  stoichiometric  flame  with  v  =  0  we  can 
also  make  use  of  Eqs.  4.9,  4.10,  and  6.21  to 
obtain  the  following  explicit  expression  for  pi¬ 
rn  terms  of  original  system  parameters 

2'^ /l,,Le^m- 

“  9=[0.6307m=  -  1.344m  +  1] 

X  (LeH^LepJ-^- 

Len.o 

-O,  ^  ®  ^ 

(6.24) 

This  represents  the  maximum  value  of  the 
burning  rate  that  is  achieved  when  the  rate  of 
radical  production  is  large.  It  is  interesting  to 
note  that  A:,,  and  k,,,  have  opposite  effects  on 
the  burning  rate.  Since  Eq.  6.24  is  obtained  in 
the  limit  /!,-*«,  there  are  many  radicals 
present  in  the  flame  structure.  The  present 
theory  predicts  that  the  flame  will  propagate 
more  rapidly  when  the  fuel  consumes  these 
radicals  at  a  faster  rate  than  they  recombine. 

When  D  is  decreased,  for  example  by  in¬ 
creasing  the  pressure,  the  flame  speed  also 
decreases  monotonically,  as  shown  in  Fig.  2. 
Eventually,  the  chain-breaking  reaction  be¬ 
comes  dominant,  the  burning  rate  decreases  to 
zero,  and  the  flame  ceases  to  propagate.  Since 
m  ~  -2  for  the  methane-air  flames  consid¬ 
ered  here,  the  correction  term  in  Eq.  6.22  is 
always  positive,  indicating  that  the  response 
curve  is  monotonic  and  no  turning  point  be¬ 
havior  is  predicted.  This  is  in  agreement  with 


previous  numerical  results  [13]  which  have 
shown  that,  in  the  absence  of  heat  loss,  the 
flame  speed  will  decrease  to  zero  monotoni¬ 
cally  and  a  distinct  extinction  state  does  not 
e^ist. 

For  the  two  limiting  cases  discussed  above 
we  are  able  to  determine  the  explicit  depen¬ 
dence  of  the  stoichiometric  burning  rate  on 
pressure.  In  particular,  when  D  -*  *  we  note 
from  Eq.  6.24  that  pv  varies  with  the  square 
root  of  pressure.  Thus  while  the  burning  rate 
increases  with  pressure,  the  burning  velocity  is 
seen  to  decrease.  For  the  other  limiting  case, 
Eq.  6.22  for  which  the  velocity  approaches 
zero,  pv  has  the  explicit  form  (with  i/  =  0) 

Pf-  “■  •  «  A  . 

""I  C^H^LCh, 0.^(1  + 


2A]Uo,  15' 


In  this  case,  both  burning  rate  and  velocity 
decrease  as  pressure  is  increased.  We  also  ob¬ 
serve  from  the  above  expression  that  an  in¬ 
crease  in  y4,,  results  in  a  decrease  in  pv.  This 
is  in  contrast  to  the  strongly  burning  limit,  Eq. 
6.24,  for  which  the  chain-propagating  reaction 
was  found  to  further  enhance  burning. 


6.2.  Structure  Equation  for  Fuel-Lean  Mixtures 

The  analysis  for  the  fuel-lean  mixture  is  nearly 
identical  to  the  stoichiometric  case  discussed 
above,  so  we  omit  most  of  the  details.  The 
form  of  the  structure  equations  remains  the 
same  as  in  Eqs.  6.12-6.13  and  6.19-6.20.  How¬ 
ever,  to  leading  order  we  now  have  (A’q^)”  “ 
Xq,  f,  “  2Leo|(<^~'  -  1).  Thus  the  parameters 
A  and  y,  given  by  Eqs.  4.9  and  6.21,  respec¬ 
tively,  are  independent  of  the  burning  rate.  An 
expression  for  D  is  obtained  by  inserting  Eq. 
5.21  into  Eq.  5.22  which  yields  the  dependence 
D  ~  {pv)~*^'.  It  follows  from  Eq.  6.21  that 
the  rate-ratio  parameter  D  also  exhibits  this 
dependence  on  burning  rate. 

In  the  limit  ^  -♦  0  we  again  obtain  Eq.  6.17 
as  the  expression  that  determines  the  burning 
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rate.  When  v  =  0  we  can  use  Eqs.  5.21  and 
5.22  to  write  the  following  explicit  expression 
for  pv  in  terms  of  all  other  system  parameters 

(5 /4)^^^(  ..4 , 1  Le 

''(2c,»',LeH,o^H.,o,*(l +  «)'''/ 

(6.26) 

The  burning  rate  is  seen  to  increase  with  pres¬ 
sure  as  pr  and  thus  burning  velocity 

decreases  with  pressure.  We  also  note  the  di¬ 
rect  dependence  on  equivalence  ratio  pi  ~ 
(4>~'  —  In  addition,  pi  is  proportional 

to  A'l^*  with  Ajj  -  exp(-£,,/R''f “),  where 
T°  decreases  with  equivalence  ratio.  Thus  the 
burning  rate  decreases  as  we  move  further 
away  from  stoichiometry  (for  1  -  =  0(c) 

the  analysis  of  the  preceding  subsection  is  ap¬ 
propriate).  Reactions  1;,,  and  k,„  are  seen  to 
have  opposite  effects  on  the  burning  rate,  in 
agreement  with  the  stoichiometric  situation  for 
this  limit.  Again  we  find  that  the  flame  propa¬ 
gates  faster  when  the  rate  of  fuel  consumption 
exceeds  the  rate  of  radical  recombination. 

In  the  limit  that  radical  production  is  very 
slow,  that  is,  ,  -*  0,  we  again  obtain  Eq.  6.22. 
For  the  present  case,  however,  D  depends 
explicitly  on  pi  so  that  the  leading  order  value 
D  =  15/8  does  provide  an  expression  for  the 
burning  rate.  In  particular,  this  equation  can 
be  written  in  the  form 

..  8(Leo^^o..^)V./^ 

AVf..Po-v/j^(LeH,)^'^  r' 

5c,H>LeH^oArH^o.»(l  +  ’ 

(6.27) 

so  that  pi  ^  A^^'  when  A,  -*  0.  Again  we 
find  the  dependence  on  pressure  to  be  pi  - 
As  before,  although  pi  -  (Xq.,/,)' 

-  1)‘,  the  Arrhenius  dependence  of  pi 
on  will  cause  the  burning  rate  to  decrease 


as  the  mbcture  is  made  more  lean,  so  long  as  ^ 
is  bounded  away  from  unity.  From  definition 
5.21  6  is  found  to  decrease  as  the  lean  limit  is 
approached,  resulting  in  a  lower  reaction  zone 
temperature  (as  detenrined  by  solving  Eqs. 
2.10  and  6.4  simultaneously)  and  consequently 
a  lower  burning  rate.  We  remark  that  the 
burning  rate  is  small  in  this  limit,  and  the 
flame  will  probably  have  been  extinguished  by 
the  presence  of  unavoidable  heat  loss  before 
the  limiting  value  6.27  can  be  reached.  This 
point  is  currently  being  investigated. 

When  I'  0  the  only  place  the  burning  rate 
appears  in  the  stnicture  equation  is  in  the 
form  or  which  is  plotted 

as  a  function  of  y  in  Fig.  3.  This  curve  can  be 
interpreted  as  showing  the  velocity  of  fuel-lean 
flames  as  a  function  of  the  ratio  of  branching 
to  propagating  reactions.  The  response  is  seen 
to  be  monotonic,  and  the  burning  rate  ap¬ 
proaches  the  limiting  values  Eqs.  6.26  and  6.27 
as  y  -♦  *  and  0,  respectively. 

7.  RESULTS  AND  DISCUSSION 

Equations  6.19  and  6.20  were  solved  numeri¬ 
cally  to  compute  the  burning  velocity  as  a 
function  of  pressure  and  equivalence  ratio  for 
both  stoichiometric  and  lean  mixtures.  Appro¬ 
priate  values  for  all  parameters  were  the  same 
as  those  used  by  Peters  and  Williams  [4],  and 
Seshadri  and  Peters  [5].  In  particular,  the  value 
of  A/Cp  ^  was  taken  to  be  A/Cp  *  2.58  x 
10*^(f/Ti,)“’  g/cm  s,  and  the  Lewis  numbers 


y 


Fig.  3,  Plot  of  representing  burning  rite  for 

fuel-lean  flames,  versus  the  parameter  y.  Curve  is  gener¬ 
ated  by  solving  equations  6.19  and  6.30. 
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were  assigned  the  values  Le^  •*  0.97,  Leo  •= 

1.1,  Le„,o  *  0.83.  Leco,  -  1.39,  Le„^  -  5.3, 
LeH  *  0.18,  and  Leco  “  Lll.  The  average 
molecular  weight  was  assumed  constant  equal 
to  27.62  kg/kmol,  and^  the  ambient  tempera¬ 
ture  was  taken  to  be  -  300  K.  The  kinetic 
data,  given  in  [4]  and  [5],  is  *  2.0  x  10’^ 
cm^/mol  •  s,  Aij  *•  2.2  x  10*  cm^/mol  •  s,  £, 
-  16800  cal/mol,  =  8750  cal/mol.  The 
reaction  rate  k,,  has  the  form  k,,  > 
.4„r^  exp( -£„/£“/)  where  £“-1.987 
cal/mol  K  (£  -  82.05  atm  cm-’ /mol  K)  is  the 
gas  constant.  Since  k,/k,,  has  a  weak  temper¬ 
ature  dependence,  we  have  approximated  this 
reaction  rate  by  k„  -  >4, i  exp( -£,/£“ f). 
Consistent  with  this  approximation,  we  have 
defined  an  effective  pre-exponential  term 

=^*ii^r^exp[(£,  -£,,)/£“7;„],  where 
we  have  chosen  1600  K.  The  remaining 
reaction  rates  and  equilibrium  constants  are 
given  by  -  1.22  cxp(3067/f),  a  - 

3.86exp(-3652/f),  £3  -  0.216  exp(7658/f), 
k,o  “  1.51  X  f'  ^  exp(758/£“f). 

The  procedure  we  used  to  determme  ''e 
burning  velocity  and  the  temperature  f  “  is  as 
follows.  For  a  given  equivalence  ratio,  the  un- 
bumed  fuel  mass  fraction  is  given  by  „  - 
(17.16  -t-  4>),  and  a  thermochemical  calcula¬ 
tion  yields  the  value  of  Equations  6.19, 
6.20,  and  6.4,  which  can  be  expressed  in  terms 
of  original  variables  as 

ft  -  f“  -  eZ(0)(f“  -  t),  (7.1) 

are  solved  simultaneously  to  determine  f  “  and 
V  as  a  function  of  p.  It  actually  proved  easier 
to  solve  for  v  and  p  in  terms  of  r“,  which  was 
first  assigned  a  value  close  to  fj,  and  then 
gradually  decreased.  For  a  fixed  value  of  f  a 
value  of  €  is  immediately  obtained  from  Eq. 

7.1.  The  solution  of  Eqs.  6.19  and  6.20  is  used 
to  determine  the  pressure  p,  and  the  velocity  is 
then  found  using  the  appropriate  definition  of 
t  given  in  Section  5. 

In  Fig.  4  we  have  plotted  v  and  f  “  as  a 
function  of  pressure  for  stoichiometric  mix¬ 
tures  assuming  partial  equilibrium  of  the  wa¬ 
ter-gas  shift  ( 1/  -  0).  At  p  -  1  atm,  our  analy¬ 
sis  yields  a  value  for  the  burning  velocity  of 
t  -  36  cm/s.  in  reasonable  agreement  with 
experimentally  measured  values  [14,  15].  The 


Fig.  4.  Velocity  and  fuel  consumption-zcne  temperature 
as  a  funaion  of  pressure  for  stoichiometric  flame  with 
partial  equilibrium  of  water-gas  shift.  Dashed  line  indi¬ 
cates  results  obtained  using  rate-ratio  asymptotics. 


velocity  .is  seen  to  decrease  as  pressure  is  in¬ 
creased,  again  in  agreement  with  observations. 
However;  the  predicted  drop  is  much  more 
rapid  than  experiments  indicate,  and  the  veloc¬ 
ity  approaches  zero  when  the  pressure  ap¬ 
proaches  a  critical  value  p^  26  atm.  As 
pointed  out  in  Ref.  4,  this  poor  agreement  at 
elevated  pressures  is  due  to  deficiencies  in  the 
kinetic  scheme.  Better  agreement  can  in  fact 
be  obtained  by  including  more  reactions  from 
the  starting  mechanisms  [4]. 

The  dashed  lines  in  Fig.  4  are  the  results 
obtained  by  setting  y  -  0  in  Eq.  6.19.  Recall 
that  in  this  limit  the  structure  problem  6.19 
and  6.20  reduces  to  the  form  found  using  rate- 
ratio  techniques,  and  we  will  refer  to  this  as 
the  rate-ratio  limit.  The  rate-ratio  approach  is 
seen  to  predict  somewhat  higher  velocities  than 
our  theory,  most  notably  at  lower  pressures.  At 
p  -  1  atm  the  rate-ratio  analysis  predicts  a 
value  for  the  burning  velocity  of  u  -  51  cm/s. 
The  slight  difference  in  this  number  and  the 
value  quoted  in  Refs.  4  and  6  is  due  to  the  way 
in  which  they  chose  to  represent  the  reaction 
rate  terms.  As  pressure  is  increased  the  dif¬ 
ference  in  the  two  theories  is  seen  to  diminish, 
and  the  curves  converge  as  p  ->  p,,  u  -*  0.  As 
discussed  earlier,  agreement  between  the  two 
theories  in  this  limit  is  not  surprising  since  the 
value  y  —  0  corresponds  to  f  -»  0  in  our  the¬ 
ory. 

Figure  5  shows  the  correction  to  the  burning 
rates  for  stoichiometric  mixtures  when 
nonequilibrium  effects  are  included.  For  both 
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Fig.  5.  Burning  rate  as  a  function  of  pressure  for  stoichio¬ 
metric  flame  illustrating  modification  due  to  nonequilib¬ 
rium  of  the  water-gas  shift.  Solid  line  indicates  results 
obtained  using  rate-ratio  asymptotics. 

theories,  nonequilibrium  effects  are  seen  to 
cause  a  reduction  in  the  burning  rate.  Again 
the  difference  in  the  curves  is  seen  to  be 
greatest  at  lower  pressures.  Partial  equilibrium 
of  reaction  II  becomes  a  better  assumption  as 
pressure  is  increased,  and  the  curves  converge 
as  p  -» p,. 

Overall,  the  agreement  between  the  rate- 
ratio  approach  and  the  present  theory  is  seen 
to  be  quite  favorable  for  premixed  methane -air 
flames.  In  both  cases  comparison  with  exp>eri- 
mental  observations  is  satisfactory,  especially 
at  lower  pressures.  The  reason  for  the  good 
agreement  between  the  two  theories  is  that  the 
parameter  y  appearing  in  our  structure  equa¬ 
tion  6.19  is  typically  quite  small.  For  the  wide 
range  of  parameters  considered  in  our  compu¬ 
tations,  rarely  did  y  approach  unity.  However, 
the  present  analysis  can  be  extended  to  other 
hydrocarbon- air  mixtures  with  the  same  basic 
structure  for  which  larger  values  of  y  are  pos¬ 
sible.  From  the  definitions  of  y  and  A,  Eqs. 
6.21  and  4.9,  respectively,  any  mixture  change 
that  results  in  a  slower  propagation  reaction 
rate  or  an  increased  O,  concentration  will  also 
yield  larger  values  of  y.  In  those  cases,  the 
structure  equations  6.19-6.20  may  yield  more 
accurate  results  than  when  rate-ratio  tech¬ 
niques  are  used  in  the  fuel  consumption  zone. 
An  additional  advantage  of  the  present  ap¬ 
proach  is  the  fact  that  we  have  retained  the 
exponential  nonlinearity  of  the  reaction  rate 
term  in  the  fuel  consumption  layer.  This  should 
enable  us  to  study  the  response  of  these  flames 
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to  small  perturbations  such  as  heat  loss  and 
stretch.  These  points  are  currently  being  inves¬ 
tigated. 

The  velocity  and  the  burning  rate  as  com¬ 
puted  from  the  present  theory  with  ^  »  1  are 
plotted  in  Fig.  6,  where  results  from  both  par¬ 
tial  equilibrium  and  nonequilibrium  .calcula¬ 
tions  are  shown.  Although  the  velocity 
decreases  with  an  increase  in  pressure,  the 
burning  rate  is  seen  to  first  increase  before 
reaching  a  maximum  value  and  then  dropping 
off.  This  trend  has  been  observed  both  numeri¬ 
cally  and  experimentally  for  dilute  mixtures 
[16],  although  preliminary  data  for  nondilute 
mixtures  does  not  conclusively  show  this  trend. 
Measurements  need  to  be  made  over  a  wider 
range  of  pressures  to  determine  whether  the 
trend  predicted  by  the  present  theory  can  be 
realized  for  non-dilute  mixtures. 

In  Fig.  7  we  compare  experimentally  mea¬ 
sured  velocities  and  burning  rates  to  values 
predicted  by  the  present  asymptotic  theory  and 
by  numerical  calculations  with  the  detailed 
C2*chain  mechanism  discussed  in  Ref.  17.  The 
experimental  measurements  were  made  in  our 
laboratory  .using  the  counterflow  twin-flame 
technique.  The  numerics  are  seen  to  agree 
quite  well  with  the  experimental  results  for  the 
entire  range  of  pressures  up  to  8  atm.  The 
asymptotic  results  show  good  agreement  up  to 
about  4  atm,  after  which  velocity  and  burning 
rate  drop  off  too  quickly.  As  mentioned  ear¬ 
lier,  better  agreement  can  be  obtained  by  in¬ 
cluding  more  elementary  reactions  from  the 
starting  mechanism. 


Fig.  6.  Velocit)'  and  burning  rate  shown  as  a  function  of 
pressure  for  stoichiometric  flame.  Solid  line  denotes  values 
obtained  when  partial  equilibrium  of  the  water-gas  shift  is 
assumed. 
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Fig.  7.  Comparison  of  experimentally  measured  velocity 
and  burning  rate  with  values  obtained  numerically  and 
with  present  asymptotic  theory.  Solid  diamonds  denote 
velocity. 


The  temperature  at  the  leading  edge  of  the 
flame  7”  and  the  burning  velocity  of  near- 
stoichiometric  flames  were  also  computed  for 
several  values  of  equivalence  ratio.  As  the  mix¬ 
ture  was  made  more  lean,  the  characteristic 
temperature  of  the  fuel  consumption  zone  f  ® 
was  found  to  decrease,  as  expected.  However, 
computed  burning  rate  and  velocity  profiles 
first  show  a  slight  increase  until  ^  0.85,  after 
which  pi'  decreases  with  further  decreases  in 
equivalence  ratio.  If  we  examine  the  explicit 
expression  for  pr,  found  by  using  Eqs.  5.10 
and  7.1, 


( pvf 


'  3(f, -f")  V 


Ayf.,p-HTLeo,LeH:)'''V^ 


V  1 

4v/2(l  +  a)’^*/’ 


(7.2) 


we  see  that  although  the  difference  between 
and  fo  decreases  with  decreasing  f°,  the  other 
quantities  in  Eq.  7.2  increase.  These  factors 
lead  to  the  increase  in  velocity'  as  <5  is  de¬ 
creased  in  the  range  (1.0,  0.85).  We  remark 
that  the  same  incorrect  trend  is  obtained  if  y 
is  set  to  zero  in  our  structure  equation  6.19, 
and  therefore  this  erroneous  prediction  is  not 
due  to  the  activation-energ\’  asymptotic  tech¬ 
niques  that  we  employed  in  the  fuel  .consump¬ 


P  (8tm) 

Fig.  8.  Velocity  and  burning  rate  of  lean  flame  as  a  func¬ 
tion  of  pressure  for  three  different  values  of  equivalence 
ratio. 

tion  zone.  Rather,  the  discrepancy  is  most  likely 
due  to  deficiencies  in  the  simplified  kinetic 
scheme.  The  increase  is  veiy  slight  (roughly 
10%),  suggesting  that  higher-order  terms  that 
we  neglected  can  reverse  this  trend.  Indeed, 
Seshadri  and  Gottgens  [6]  and  Bui-Pham  et  al. 
[7]  have  included  many  more  intermediate  re¬ 
actions  in  their  rate-ratio  asymptotic  analyses 
of  lean  methane -air  flames,  and  they  report 
good  agreement  with  experiments,  cf.  Fig.  13 
of  Ref.  6.  The  inclusion  of  these  additional 
reactions  in  our  analysis  will  also  yield  the 
appropriate  trend,  although  the  additional  pa¬ 
rameters  that  enter  into  the  analysis  greatly 
complicate  the  details.  It  is  clear  that,  if  one 
desires  to  compute  accurate  values  for  burning 
velocities  over  a  wide  range  of  parameters, 
many  intermediates  need  to  be  retained  in  the 
analysis. 

As  the  mixture  is  made  more  lean  such  that 
1  -  >  0(1),  an  abundance  of  oxygen  passes 

through  the  flame  unreacted  and  the  structure 
equation  discussed  in  Section  6.2  becomes  ap¬ 
propriate.  As  mentioned  earlier,  the  system  of 
equation!)  6.19  and  6.20  remain  the  same  but 
the  definitions  of  D  and  y  change.  In  Fig.  7  we 
show  the  velocity  of  these  lean  flames  as  a 
function  of  pressure  as  computed  for  several 
values  of  equivalence  ratio.  Now  the  flame  is 
seen  to  bum  weakest  as  the  lean  limit  is  api- 
proached  and,  near  this  limit,  the  velocity  in¬ 
creases  with  As  is  increased  further  such 
that  1  -  ^  «  0(e),  the  analysis  of  Section  5.2 
is  no  longer  valid,  the  structure  of  the  oxida¬ 
tion  layer  changes,  and  the  analysis  of  Section 
5.1  is  appropriate. 
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8.  SUMMARY  AND  CONCLUSIONS 

The  structure  of  premixed  methane-air  flames 
with  a  detailed  four-step  kinetic  mechanism 
has  been  studied  theoretically.  The  basic  flame 
structure,  first  identified  by  Peters  and  Williams' 

14] ,  remains  the  same  as  in  previous  rate-ratio 
asymptotic  analyses.  However,  in  the  present 
work  the  effective  activation  energies  of  the 
intermediate  reactions  describing  the  competi¬ 
tion  of  fuel  and  oxidant  for  radicals  are  as¬ 
sumed  to  be  sufficiently  large  such  that  activa¬ 
tion  energy  asymptotics  can  be  employed  in 
the  fuel  consumption  zone. 

Our  results  include  the  evaluation  of  the 
burning  rate  in  terms  of  pressure,  equivalence 
ratio,  and  the  ratio  of  the  competing  rates  of 
the  branching  reaction  to  the  termination  reac¬ 
tion.  Although  limitations  of  the  simplified  ki¬ 
netic  scheme  were  noted  in  the  analysis  of 
near-stoichiometric  flames,  for  the  most  part 
the  trends  are  found  to  be  in  qualitative  agree¬ 
ment  with  experimental  observations,  cf.  (14, 

15] .  For  the  adiabatic  conditions  considered 
here,  the  burning  velocity  does  not  exhibit 
turning-point  behavior  (i.e.,  extinction),  in 
agreement  with  numerical  computations  [13]. 

We  remark  that  although  the  activation  en¬ 
ergies  of  the  intermediate  reactions  are  not 
significantly  large,  our  results  suggest  that  they 
are  sufficiently  large  to  permit  the  use  of  acti¬ 
vation  energy  asymptotics.  Our  results  agree 
quite  well  with  those  obtained  using  rate-ratio 
asymptotics,  and  in  fact,  the  structure  equation 
found  with  rate-ratio  techniques  can  be  recov¬ 
ered  as  a  limiting  case  of  the  equation  that  we 
have  derived.  Specifically,  by  setting  y  «  0  in 
£q.  6.19,  the  equation  first  derived  in  Ref.  4  is 
obtained.  For  the  premixed  methane -air 
flames  of  primary  interest  here,  the  value  of  y 
is  typically  quite  small.  However,  it  is  possible 
to  consider  other  hydrocarbons  or  to  modify 
the  mixture  in  such  a  way  as  to  obtain  larger 
values  of  y  without  changing  the  basic  flame 
structure.  For  such  situations  there  will  cer¬ 
tainly  be  less  quantitative  agreement  between 
the  two  theories.  As  discussed  in  Section  6, 
qualitative  differences  will  also  arise  for  larger 
values  of  y.  For  example,  the  RRA  theory 
predicts  that  Ic,,  is  always  terminating,  while 
the  present  theory  indicates  that  its  role  is 
reversed  for  larger  values  of  y.  One  advantage 


of  the  method  proposed  here  is  that  the  non¬ 
linear  reaction  rate  term  is  retained  in  the  fuel 
consumption  zone.  Experience  with  one-step 
chemistiy  suggests  that  such  nonlinearities  can 
realistically  capture  a  variety  of  flame  re¬ 
sponses  to  perturbations.  Although  Damkdhler 
number  asymptotics  and  rate-ratio  techniques 
are  needed  to  analyze  many  of  the  layers  within 
these  complicated  flame  structures,  by  employ¬ 
ing  activation  energy  asymptotics  in  the  fuel 
consumption  layer,  perhaps  we  can  better  study 
the  response  of  these  hydrocarbon  flames  to 
heat  loss,  stretch,  etc. 

To  illustrate  our  method  we  have  used  the 
most  basic  kinetic  mechanism,  which  repre¬ 
sents  a  significant  simplification  from  the  origi¬ 
nal  starting  mechanism,  and  consequently  there 
is  a  certain  degree  of  error  inherently  present. 
However,  more  intermediate  reactions  can  be 
retained,  and  in  addition  several  assumptions 
regarding  the  relative  thickness  of  the  various 
zones  can  be  relaxed  in  order  to  gain  increased 
accuracy.  Seshadri  (cf.  Refs.  5-7)  has  done  a 
great  deal  of  work  in  refining  this  mechanism 
to  obtain  much-improved  results.  The  results 
of  those  studies  agree  favorably  with  the  de¬ 
tailed  numerical  calculations  of  Smooke  and 
Giovangigli  [2]  and  Egolfopoulos  et  al.  [15],  as 
well  as  the  experimental  measurements  re¬ 
ported  in  Ref.  15.  The  same  improvements  can 
be  made  in  the  present  study  to  obtain  better 
agreement  with  numerics  and  experiments. 

As  a  final  comment,  we  remark  that  the  four 
global  reactions  used  in  the  present  study  can 
also  be  used  to  investigate  the  structure  of 
methane-air  diffusion  flames,  as  has  been  done 
using  rate-ratio  asymptotics  (cf  Refs.  18-20). 
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Analysis  of  Thermal  Ignition  in  the  Supersonic  Mixing  Layer 
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Ignition  in  a  buninar  tuptnonic  mixing  layer  between  two  parallel  etreams  of  initially  teparated  reactants  b 
studied  both  numerically  and  through  the  use  of  large  activation  energy  asymptotics.  The  asymptotic  analysb 
provides  a  description  of  ignition  characteristics  over  the  entire  range  of  system  parameters.  In  partfeular,  k  b 
demonstrated  that,  for  small  values  of  viscous  heating,  the  ignitioo  distance  acalm  approximately  linearly  with 
the  freestream  Mach  number,  whereas  for  large  viscous  beating  h  decreases  rapidly  due  to  the  tcmperatifre-scn- 
sitive  nature  of  the  reaction  rate.  Tbb  indicates  the  potential  of  using  local  flow  retardation  to  enhance  ignitioo 
rather  than  relying  solely  on  external  heating.  The  asymptotic  analysb  further  Mentifles  several  distinct  ignition 
situations,  yielding  resuhs  that  compare  well  with  those  obtained  flum  the  full  numerical  cakulatioa.  The  effeeb 
of  flow  nonsimllarity  arc  abo  assessed  and  are  found  to  be  more  prominent  for  the  mixing  layer  flow  in  compari¬ 
son  to  the  flat-plate  ccmflguration  studied  previously. 


Introduction 

HE  study  of  combustion  within  supersonic  boundary-layer 
flows  has  generated  much  recent  interest  due  to  its  relevance 
in  the  development  of  scramjet  engines.  The  most  essential  feature 
of  supersonic  boundary-layer  combustion,  in  contrast  to  the  sub¬ 
sonic  situation,  is  that  the  high-speed  flow  contains  a  considerable 
amount  of  kinetic  energy  that  can  be  convened  to  thermal  energy 
through  viscous  dissipation,  thus  facilitating  ignition  of  the  com¬ 
bustible  mixture. 

In  a  recent  snidy'  we  investigated  the  effect  of  viscous  heating 
on  the  ignition  of  a  supersonic  stream  of  premixed  combustible 
over  a  flat  plate,  using  full  numerical  compuution  and  large  acti¬ 
vation  energy  asymptotics.  Several  distinct  ignition  regimes  were 
identified,  d^nding  on  the  relative  magnitude  of  the  external  heat 
transfer  from  the  hot  boundary  and  the  internal  heat  generation  due 
to  viscous  heating.  The  analysis  of  each  regime  was  shown  to 
properly  capture  the  ignition  characteristics. 

Whereas  the  flat-plate  configuration  studied  previously*  pro-  • 
vides  a  simple  model  problem,  many  qrplications  of  interest,  most 
notably  the  scramjet  engine,  involve  the  compressible  mixing  layer 
between  fuel  and  oxidizer  supplies.^  In  contrast  to  the  flat  plate, 
the  problem  is  nonpremixed  in  nature;  therefore  mixing  of  the  re¬ 
actants  must  occur  to  achieve  ignition.  Furthermore,  the  internal 
generation  of  heat  in  the  presence  of  a  shear  layer  will  vary  de¬ 
pending  on  the  difference  in  velocities  of  the  two  parallel  streams. 
Thus  the  ignition  characteristics  of  the  mixing  layer  are  expected 
to  be  different  from  that  of  the  supersonic  flat-plate  flow,  for 
which  the  no-slip  condition  at  the  wall  always  resulu  in  large 
amounts  of  viscous  heating.  Therefore,  in  this  paper  we  extend  our 
previous  work*  and  consider  a  mixing  layer  configuration  to  pro¬ 
vide  a  comprehensive  treatment  on  the  subject  of  ignition  within 
supersonic  boundary-layer  flows. 

The  effect  of  viscous  heating  on  the  ignition  in  a  supersonic 
mixing  layer  was  first  considered  asymptotically  by  Jackson  and 
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Hussaini,^  in  which  nearly  equal  ffeestxeam  temperatures  and  ve¬ 
locities  were  assumed.  Ju  and  Niioka*  allowed  for  an  0(1)  temper¬ 
ature  difference  between  the  two  fteestreams,  nonetheless  restrict¬ 
ing  the  analysis  to  small  velocity  differences.  A  pseudomultittep 
chemistry  model  was  also  used.  A  more  thorough  treatment  of  0(1) 
viscous  heating  in  the  mixing  layer  was  subsequently  made  1^ 
Grosch  and  Jackson,^  in  which  Jackson's  and  Hussaini’s  psevious 
aiudysis^  was  extended  to  tender  it  valid  over  a  wider  range  of  pa¬ 
rameters  Characterizing  the  strength  of  the  ignition  source.  A  nu¬ 
merical  and  asymptotic  analysis  of  the  evolution  from  a  weakly  re¬ 
active  state  to  the  fully  developed  laminar  diffusion  flame  was  also 
provided,  yielding  a  more  accurate  description  of  die  diffusion 
flame  structure.  It  was  further  emphasized  dut,  because  combus¬ 
tion  in  the  mixing  layer  can  introduce  additional  flow  instabilities,* 
which  in  turn  enhaiKes  reactant  mixing,  it  is  of  fundamental  im¬ 
portance  to  investigate  the  effect  of  chemical  reaction  on  the  flow- 
field.  Furthermore,  a  good  understanding  of  ignition  and  flame 
propagation  in  an  idealized  laminar  mixing  layer  provides  useful 
information  for  more  complicated  studies  involving  chemicaUy  re¬ 
acting  flows. 

The  present  study  extends  the  work  of  Grosch  and  Jackson*  by 
treating  all  possible  weakly  reactive  regimes  that  can  exist  in  su¬ 
personic  mixing  layers  prior  to  the  point  of  thermal  rutuway.  This 
includes  the  case,  not  considered  in  Ref.  5,  in  which  the  heat  for  ig¬ 
nition  generated  by  the  hot  stream  and  through  viscous  heating  are 
comparable.  Furthermore,  for  each  ignition  situation,  we  derive 
explicit  expressions  for  the  minimum  ignition  distance,  which  is 
then  mapp^  out  over  the  complete  range  of  system  parameters. 
The  results  compare  well  with  numerical  computations.  The  sig¬ 
nificance  of  flow  nonsimilaiity  is  discussed,  and  comparisons  with 
the  flate-plate  case*  are  also  tnade. 

In  the  following  section,  the  model  problem  is  formulated  and 
the  distinct  ignition  situations  ate  classified.  Next,  die  numerical 
solutions  for  the  full  system  of  equations  are  presented.  We  then 
proceed  with  the  asymptotic  analysis,  and  the  results  are  compared 
with  numerical  pr^ictions.  Finally,  we  add  further  discussions 
and  concluding  remarks. 

Formulation 

As  shown  in  Fig.  1,  we  consider  a  stream  of  reactant  1  with  re¬ 
actant  mass  fraction  f , density  velocity  and  teiiqieta- 
ture  t„,  flowing  parallel  to  a  stream  of  reactant  2,  with  f  2.^, 

and  f...  We  assume  a  zero  pressure  gradient  existt  across 
die  boundary  layer  such  that  the  inert  solution  has  a  self-similar 
form.  The  two  initially  separated  reactants  diffuse  into  one  anodi- 
er,  thus  permining  ignition  to  occur  when  the  temperature  is  ele¬ 
vated  to  a  sufficiently  high  value.  Chemical  reaction  it  weak  near 
the  leading  edge  of  the  mixing  layer,  and  it  evolves  gradually  in 
the  downstream  direction.  Eventu^ly  an  abrupt  temperature  rise. 
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coiiesponding  to  thermal  runaway,  occurs  at  a  certain  location 
within  the  mixing  layer.  We  are  interested  in  analyzing  the  reac¬ 
tion  ztme  structure  up  to  the  point  of  ignition  to  obtain  explicit  ex¬ 
pressions  for  the  minimum  ignition  distance  as  a  function  of  all  the 
system  parameters.  We  assume  that  a  one-step,  overall  reaction, 
with  an  Arrhenius  temperature  dependence,  takes  place  between 
species  1  and  2  to  form  the  product.  We  further  assume  constant 
properties  such  as  specific  heat  c^.  density-weighted  viscosity  p  p , 
and  unity  values  for  the  Lewis  and  Prandtl  numbers.  The  conserva¬ 
tion  equations  for  momentum,  mass  fraction  of  species  i.  and  static 
energy  in  the  supersonic  mixing  layer  are  then  given  by^ 

/“•  +  //"=0  (1) 

/IV  /IV 

3^  =  C>',rT'-'’-‘'*‘exp(-Vn  (2) 

-2hI/''(t1)]V(1-X)^  (3) 


where  (  )'  denotes  d/dq,  and  the  preceding  equations  are  to  be 
supplemented  by  the  ideal-gas  equation  of  state.  The  nondimen- 
sional  temperature  T  and  the  mass  fraction  for  species  /,  y,  are  de¬ 
fined  in  terms  of  the  original  variables  by  T  =  f  2._  and  T,  = 
f’i  /O/^  j.—  respectively.  Here,  g  is  the  heat  of  reaction  per  unit 
mass  of  species  2,  and  o,  is  the  stoichiometric  mass  ratio  of  each 
species  relative  to  species  2.  such  that  Oj  =  1 .  Other  parameters  in¬ 
clude  the  nondimensional  activation  temperature  =  CpiEJF^)! 

j.^,  ard  the  viscous  heating  parameter  p  ■  (1/2)(y  - 
(1  -  where  y^CpICy,  1  =  and 


Af-  *  nyJyPtlTw 


is  the  Mach  number  of  the  upper  freestream.  The  variables  ^  and  r\ 
■re  the  normalized  Howarth-Dorodnitsyn  variables  defined  as 


■nd/is  the  stream  function  v(x,  v)  normalized  as 
/(Tl)-V(jr,y)/(2^.p.C^)''= 


(4) 


(5) 


(6) 


Here  p  and  q  represent  the  reaction  orders  of  species  1  Kid  2,  re¬ 
spectively,  and  r  is  the  temperature  exponent  of  the  pre-exponen¬ 
tial  reaction  rate  term.  In  Eq.  (4).  t|(,  denotes  a  shift  of  the  origin;  it 
may  remain  indeterminate,'  but  we  are  free  to  translate  the  coordi¬ 


nate  such  that  q  s  0  at  the  origin.  We  note  that  the  streamwise  co¬ 
ordinate  ^  is  proportional  to  the  frequency  factor  B  and  hence  can 
be  interpreted  as  a  reduced  Damkfihler  number,  representing  the 
ratio  of  a  characteristic  flow  time  to  a  characteristic  reaction  time. 
In  these  expressions,  H'  is  the  molecular  weight  of  species  i,  W  the 
average  molecular  weight,  and  x  and  y  the  i^ysical  coordinates 
parallel  and  normal  to  the  stream.  The  last  term  in  the  energy  equa¬ 
tion  (3)  represents  the  contribution  due  to  viscous  heating,  i.e.,  die 
amount  of  kinetic  energy  that  is  converted  to  thermal  energy  with¬ 
in  the  mixing  layer. 

Equations  (1-3)  are  u>  be  solved  subject  to  the  boundary  condi¬ 
tions 


/'(-)=  1.  /(0)  =  0.  /'(— )  =  X  (7) 

yi(«.;)=yi,~.  y,(— ,;)=o  (S) 

yj(~.;)=o,  y2(-«,,;)  =  i  (9) 

r(«-.;)=T-.  r(-«,;)  =  r-.-r--P  (lO) 


It  may  be  noted  that  Eqs.  (1)  and  (7),  which  govern  the  flowfield, 
are  decoupled  from  the  energy  and  species  equations  and  can  be 
solved  independently  to  yield  the  self-similar  Blasius  profile  for 
mixing  layer  flows.’  ’^  Without  loss  of  generality,  we  consider  0  i 
X  ^  1,  so  that  the  upper  stream  is  always  at  a  higher  velocity.  The 
case  X  =  1,  corresponding  to  a  uniform  parallel  flow  of  initially 
separated  reactants  with  no  viscous  heating,  has  been  studied 
LiAdn  and  Crespo."  In  the  main  text,  we  consider  the  upper  stream 
to  be  the  honer  of  the  two,  i.e.,  P  >  0.  For  the  case  P  <  0,  ignition 
can  occur  near  the  slower  stream  when  ipi2  p.  The  solution  for 
this  case  can  be  obuined  by  redefming  tte  coordiiute  system,  as 
discussed  in  Appendix  A  for  completeness. 

Coupling  functions  that  relate  the  mass  fractions  of  species  1 
and  2  to  temperature  are  found  by  integrating  aRnopriaic  chemis¬ 
try-free  linear  combinations  of  Eqs.  (2)  and  (3),  yielding 

y,  =  r.  4-  y,..  -(P + y,.-)4 + -  4)  -  r  (ii) 

yj  =  r.-(P- D^-i-p^d -4)-r  (I2) 

where  we  have  made  the  convection-free  coordinate  transforma¬ 
tion 


f  (1 -/’)/(! -X)  forXrtl 

4=  ^  0S4S1 

I  erfc  m/j2)  n  for  X  =  1 

(13) 

The  relations  (11)  and  (12),  together  with  the  previously  deter¬ 
mined  flowfield,  can  now  be  inserted  into  Eq.  (3)  to  reduce  the 
problem  to  a  single  equation  for  T. 

In  the  weakly  reactive  state,  the  reaction  term  in  Eq.  (3)  is  small 
and  can  be  neglected  in  a  first  approximation.  The  frozen  tempera¬ 
ture  7)  is  governed  by 


(14) 


which,  when  solved  subject  to  the  boundary  ccmditions  (10),  yields 

7>=T.-a5-M5^  (15) 

where  a  ■  P  -  p.  The  frozen  solution  (IS)  serves  as  the  initial 
condition  for  the  numerical  solution  and  as  the  basic  solution  for 
the  asymptotic  analysis.  As  will  be  discussed  later  in  the  asymptot¬ 
ic  analysis,  the  point  of  the  maximum  frozen  temperature  identifies 
the  ignition  location  for  systems  with  large  activation  energy  and 
varies  with  the  system  parameter  o.  which  indicates  the  relative 


343 


IM.  CHAO.  BeCHTOLO.  AND  LAW.  tfSKMALIONtTION  IN  MOONO  LAYER 


inqwftaBce  ei  either  die  external  ignition  source  lepresemed  by  3 
or  the  uitemal  source  represented  by  (1. 

In  the  next  section,  we  shall  study  the  ignition  dynamics  for  the 
preceding  t«ra  cases  through  direct  numerical  integration  of  the  en¬ 
ergy  equation  (3)  with  the  bound^  conditions  (10).  We  shall  then 
proceed  to  an  asymptmic  analysis  and  derive  explicit  expressions 
fOT  the  ignition  distance  as  functions  of  all  relevant  parameters. 

Numerical  Solutions 

Equation's  (3)  attd  (10)  are  solved  numerically  using  a  second- 
mler  finite-difference  scheme  with  implicit  determination  in  the  r) 
direction  and  marching  along  the  streamwise  coordinate  (.  The 
step  size  of  is  readjusted  as  the  solution  approaches  the  ignition 
point  to  capture  the  abrupt  temperature  rise.  In  the  present  calcula¬ 
tions  we  have  set  p  s  ^  s  r  s  1  for  simplicity.  The  parameter  values 
used  are  $  c  1.0  X  l(r  kcal/kg,  Cp  -  0.25  kcal/kg-K,  £.,  =  60  kcal/ 
mole,  and  y=  1.4,  which  are  typical  numbers  for  hydrocarbon/air 
mixtures.  Boundary  conditions  for  species  concentrations  are 
taken  to  be  f  =  0.3,  f  i,»  *  0.3  with  Oi  =  3.  which  correspond 
to  the  scaled  values  i ,  y,  _  =  0.2. 

The  streamwise  evolution  of  the  temperature  profile  within  the 
mixing  layer  is  illustrated  in  Figs.  2a  and  2b.  Figure  2a  corre¬ 
sponds  to  the  subsonic  case  in  the  absence  of  viscous  heating  (M_ 
»  0),  with  =  1000  K,  -  500  K,  and  X  =  0.  This  particular 
case  is  similar  to  the  problem  studied  in  several  previous  hot 
boundary  ignition  analyses.'^'"  It  is  seen  that  as  ^  increases,  the 
maximum  temperature  at  first  rises  gradually,  but  then  exhibits  s 
dramatic  increase,  resulting  in  the  thermal  runaway  behavior  at  ^ 
=  1.64  X  10'^.  Solutions  for  the  other  extreme  case,  for  which  vis¬ 
cous  heating  is  the  only  heat  source,  ate  shown  in  Fig.  2b.  where 
we  have  set  f «  =  f  __  *  500  K,  X  =  0,  and  Af_  =  4.472  so  that  the 


Fig.  2  Evoiuthm  of  the  temperature  profiles  as  computed  numeri- 
call)  tor  a)  the  hol-tircam  case  with  f .  >  1000  K.  f ■  SOO  K.  Af.  ■  0 
tmd  b)  the  vtacous  heating  case  srith  f .  «  f  ^  >  500  K,  M.  ■  4.471 


a) 


c 

b) 

Fig.  3  Streamwise  variation  of  the  maximum  temperature  for  a)  the 
hot-stream  case  with  same  parameter  values  as  in  Fig.  2a:  vertical  lines 
represent  asymptotic  solution,  Eq.  (24),  with  A,  leading-order  result 
only  and  B,  additional  subdominant  term  retained;  b)  the  viscoas  heat¬ 
ing  case  with  same  parameter  values  as  in  Fig.  2b:  vertical  lines  repre¬ 
sent  asymptotic  results  for  C,  fiilly  nonsimilar  solution,  Eq.  (33);  and 
D,  locally  similar  solution,  using  Eq.  (34). 


initial  temperature  profile  becomes  symmetric  about  ^  w  1/2  and 
the  maximum  temperature  takes  the  same  value  as  that  in  Fig.  2a. 
As  a  result,  the  most  rapid  temperature  rise  is  seen  to  take  place 
near  ^  w  1/2,  and  thermal  runaway  occurs  at  (  »  ].23  x  10'^. 

The  streamwise  variation  of  the  maximum  temperature  for  the 
above  two  cases  is  plotted  in  Figs.  3a  and  3b.  respectively.  The  ex¬ 
istence  of  a  critical  ignition  distance  is  cleariy  illustrated.  Figure 
3a  represents  solutions  in  the  absence  of  viscous  heating,  and 
curves  are  drawn  for  two  different  values  of  X.  with  all  other  pa¬ 
rameters  taken  to  be  the  same  as  in  Fig.  2a.  The  two  curves  demm- 
strate  that,  for  a  given  strength  of  the  ignition  source,  i.e.,  for  a 
fixed  value  of  p,  increases  as  the  flow  velocity  of  the  slower 
stream  increases,  'niis  indicates  that  the  ignition  distance  is  pro¬ 
portional  to  the  velocity  at  the  ignition  location.  It  should  be  point¬ 
ed  out.  however,  that  the  dependence  of  on  X  is  quite  weak.  This 
is  physically  reasonable  given  that  ignition  occurs  near  the  hot 
boundary,  4  =  0,  which  moves  with  velocity  fi.  (see  Ref.  10). 
Thus  an  increase  in  £_  results  in  only  a  slight  increase  in  the  ve¬ 
locity  field  near  the  opposite  boundary  at  4  =  0.  The  asymptotic 
analysis  of  the  next  sectitm  also  reveals  the  weak  dependence  of 
on  X  for  this  particular  ignition  situatitm. 

To  assess  the  importance  of  reuining  nonsimilar  effects,  we 
have  also  solved  die  problem  by  neglecting  the  d/d4  term  in 
Eq.  (3).  The  results  are  ploned  in  Fig.  3a.  where  the  curves  termi¬ 
nate  at  a  point  where  they  approach  an  infinite  slope.  Thus  for  lo¬ 
cally  similar  flows,  ignition  is  identified  as  the  turning  point  trf 
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these  curves,  in  contrast  to  the  full  evolutionary  problem  that  ex¬ 
hibits  thermal  runaway.  These  locally  similar  solutions  are  seen  to 
underestimate  the  ignition  distance,  consistent  with  previous  re¬ 
sults.*  The  present  results,  however,  show  a  more  significant  devi¬ 
ation  between  the  actual  value  of  and  that  obtained  with  the 
local-similarity  approximation.  This  can  be  explained  by  compar¬ 
ing  the  relative  magnitudes  of  the  diffusion  and  convection  terms 
for  both  flow  configurations.  In  the  mixing  layer,  the  reaction  zone 
is  broader  (in  physical  scale)  than  it  is  for  the  flat  plate,  thereby  re¬ 
sulting  in  smaller  transverse  temperature  gradients  whereas  the 
stieamwise  convection  term  is  larger  due  to  the  greater  flow  veloc¬ 
ity  at  the  ignition  location.  We  note  that,  in  the  flat-plate  configu- 
ratitm,'  the  local  streamwise  velocity  at  the  ignition  point  was 
found  to  be  very  small  when  ignition  occurs  near  the  hot  plate.  A 
more  detailed  discussion  of  nonsimilarity  effects  will  be  given  in 
the  next  section. 

Figure  3b  shows  a  similar  plot  of  the  maximum  temperature 
evolution  for  the  case  of  ignition  induced  purely  by  viscous  heat¬ 
ing,  where  parameter  values  are  the  same  as  those  used  in  Fig.  2b 
When  viscous  heating  is  included,  it  is  difficult  to  compare  the  de¬ 
pendence  of  ^  on  X  in  a  straightforward  manner  as  is  done  in  Fig. 
3a,  since  variations  in  X  affect  the  viscous  heating  term,  which  in 
turn  shifts  the  maximum  temperature  location.  The  curves  shown 
in  Fig.  3b  are  for  the  case  X  *  0  specifically,  and  it  is  clear  that  the 
local-similarity  approximation  substantially  underestimates  the  ac¬ 
tual  ignition  distance. 

Asymptotic  Analysis 

In  this  section  we  study  ^s.  (1-3)  and  (7-10)  using  large  acti¬ 
vation  energy  asymptotics.  We  consider  the  weakly  reactive  sute 
and  identify  the  ignition  distance  ^  as  the  point  at  which  we  ob¬ 
serve  either  thermal  runaway  or  turning  point  behavior. 

As  mentioned  earlier,  the  flowfield  is  known.*  '®  Thus  when 
^s.  (II)  and  ( 1 2)  are  used  we  only  have  to  consider  a  single  equa¬ 
tion  for  T,  with  appropriate  boundary  conditions.  These  can  be  ex¬ 
pressed  in  terms  of  the  ^  coordinate  as 


.  -;i'ryjr'*'’-**'exp(-r,/r)  (i6) 

mo-r.,  r(i,;)«r_.  r(^.  0)  =  r/(^)  (i?) 

Here  the  frozen  profile  Tf.  given  by  Eq.  (15),  provides  the  initial 
condition  at  the  leading  edge  ((  =  0)  of  the  mixing  layer.  We  use  c 
«  r,/r,  as  the  small  parameter  in  the  analysis,  where  T,  is  the 
maximum  value  of  the  frozen  temperature. 

The  parabolic  fonn  of  the  frozen  profile  (13)  and  the  numerical 
results  shown  in  Fig.  2  suggest  that  three  distinct  regimes  of  igni¬ 
tion  should  be  considered  depending  on  the  range  of  the  parameter 
a  as  shown  schematically  in  Fig.  4.  When  a  >  0  (Fig.  4a),  the  fro¬ 
zen  temperanire  has  a  maximum  value  of  T,  ■  r«  at  %  ■  0.  This  is 
called  the  hot-stream  case,  since  the  dominant  ignition  source  is 
die  hot  temperature  of  the  freestream.  Here  the  linear  behavior  of 
the  frozen  temperature  near  4  »  0  implies  that  reaction  is  confined 
within  the  layer  of  Cfe)  thickness  provided  the  activation  energy  is 
sufficiently  large.  On  the  other  hand,  when  a  <  0  (Fig.  4c),  the 
maximum  temperature  T,  ■  T.  ♦  aV4p  is  attained  in  the  interior  of 
the  mixing  layer  at  4  ■  (I  -  p/p)/2  >  0.  In  this  viscous  heating 
case,  ignition  is  characterized  by  the  internal  viscous  heating  rep¬ 
resented  by  the  parameter  p.  The  strucnire  of  the  mixing  layer  con- 
sisu  of  two  frozen  outer  zones  of  0(  I )  extent.  0<4<4fand4f<4 
<  1.  separated  by  a  thin  reaction  zone.  The  parabolic  form  of  T. 
near  its  maximum  suggests  that  the  reaction  zone  has  0(  Ve  )  thick¬ 
ness.  Finally,  as  a  approaches  zero  (Fig.  4b),  i.e.,  when  P  and  p  are 
comparable,  neither  of  the  analyses  is  valid  and  thus  a  different  an¬ 
alytical  treatment  is  required.  In  this  intermediate  case,  the  maxi¬ 
mum  frozen  temperanire  occurs  at  4  «  0,  but  now  the  reaction  zone 


is  of  CKVr)  near  its  maximum,  which  is  broader  than  die  case 
shown  in  Fig.  4a.  We  now  proceed  to  analyze  Eqs.  (16)  and  (17) 
for  these  dir^  different  cases  in  sequence. 

Hot  Stream  Case:  a  >  0 

As  observed  in  Fig.  4a.  in  this  case  the  reaction  is  confined  to 
the  thin  layer  of  C(e)  near  4  -  0  at  which  the  frozen  temperature  at¬ 
tains  the  maximum  value  of  T..  To  study  the  structure  of  this  lay¬ 
er,  we  introduce  the  stretched  coordinate  x  oi4A  znd  seek  an 
inner  solution  of  the  form 

7',„(X.i;)*7>(X)  +  eeo(X.;)  +  CKc2)  08) 

After  substituting  diei  e  expansions  into  Eq.  (16),  and  using  the  as¬ 
ymptotic  behavior  of  the  tiowfield  (see  App^ix  B) 

(19) 

we  can  write  the  leading  order  structure  equation  as 
d^6 

*  -A  (X  -  afro)  *exp(eo  -  x)  (20) 


where 

A  =  - exp(-T,/TJ  (21) 

Equation  (20)  is  to  be  solved  subject  to  the  boundary  conditions 


C7va 

eo(0.4)*o,  ■g^(«.4)*o  (22) 

where  the  laner  condition  is  obtained  by  matching  to  the  outer  so¬ 
lution,  which  is  expanded  as  4)  =  ^(^)  ^o(^'  0  ♦  ®(6*). 

The  previous  system  is  identical  to  that  found  by  LilUui  a^ 
Crespo"  for  a  uniform  parallel  flow  of  initially  separated  reac¬ 
tants.  This  implies  that,  when  the  hot  boundary  is  the  dominant  en- 


Fig.4  Schematic  of  the  frozen  temperature  profile  in  4  coordinate  for 
a)  the  hol-strcam  case,  b)  the  intermediate  case,  and  c)  the  vtacous 
heating  case. 
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Fif.  5  Lo(>log  plot  of  the  predicted  Ignition  distance  as  a  ftinction  of 
the  frecstrcam  Mach  number  for  the  hot>stream  case. 


ergy  source,  the  ignition  characteristics  are  minimally  affected  by 
the  additional  effects  of  viscous  heating  in  the  interior  of  the 
boundary  layer.  The  only  modification  here  is  the  parameter  a, 
which  represents  the  combined  effects  of  the  temperature  differ¬ 
ence  of  the  two  streams  ^  and  the  viscous  heating  effect 
Equations  (20)  and  (22)  are  solved  numerically,  and  the  solution 
for  6o  is  found  to  have  a  turning  point  at  a  critical  Damkbhler  num¬ 
ber  A;  beyond  which  no  solution  exists.  The  value  of  A;  can  be 
evaluated  numerically  as  a  function  of  the  parameter  o,  ifter 
which  the  minimum  ignition  distance  can  be  obtained  from  the 
definition  of  A  in  Eq.  (21 ),  yielding 

Cy  =  A,  (a)  (23) 

For  q  =  1 .  Liftin  and  Crespo’ '  provided  a  useful  correlation  for  the 
ignition  Damkohler  number  as  a  function  of  the  parameter  a  in  the 
form 

A,(a)  =  2e-2(2-a)/(l -0)2  (24) 


Therefore,  the  nonsimilar  terms  can  be  expected  to  play  a  some¬ 
what  larger  role  in  the  mixing  layer  situation,  although  still  re¬ 
main  as  a  higher  order  effect. 

In  Fig.  3a  the  vertical  dotted  lines  A  and  B  represent  the  ignition 
distance  for  3/.  s  0  obtained  from  Eq.  (23).  the  line  denoted 
by  A,  only  the  leading-order  term  describing  the  behavior  of  the 
flowfield  near  the  hot  boundary  is  retained  (see  Appendix  B), 
whereas  for  line  B,  a  constant  subdominant  term  is  retained  as 
well,  as  is  done  by  LiiUbi  and  Crespo."  The  latter  case  requires 
that  U£  be  replaced  by  iJdlntJa)  in  Eq.  (23).  Although  both  re- 
sulu  are  asymptotically  equivalent,  the  resulting  values  for  the  ig¬ 
nition  distances  are  seen  to  differ  by  a  considerable  amount  This  is 
due  to  the  fact  that  the  value  of  a  used  in  the  calculation  is  quite 
small. 

The  predicted  ignition  distance  for  dre  hot-stream  case  as  a 
fimction  of  the  freestream  Mach  number  is  shown  in  Fig.  S  for  sev¬ 
eral  values  of  X.  We  have  chosen  to  plot  the  quantity  Mj^,  since  it 
is  directly  proportional  to  the  actual  ignition  distaiKe  in  physical 
length.  The  numerical  and  asymptotic  results  agree  quite  well,  ex¬ 
cept  for  large  when  viscous  heating  becomes  impoitmL  When 
Af.  is  sufficiently  large,  the  ignition  point  shifts  to  the  interior  of 
the  boundary  layer.  The  present  analysis  is  therefore  no  longer  val¬ 
id,  and  the  analysis  of  the  next  section  becomes  appropriate. 

As  a  fmal  comment  for  the  hot-stieam  case,  in  a  uniform  parallel 
flowfield,  i.e.,  X  =  1,  an  increase  in  M.  does  not  generate  addition¬ 
al  heat,  and  the  ignition  distance  simply  increases  lineariy,  as 
shown  in  Fig.  S.  In  the  presence  of  shear  between  the  two  streams, 
i.e.,  when  X  *  1,  is  seen  to  achieve  a  maximum  value,  after 
which  viscous  heating  causes  it  to  decrease,  thereby  enhancing  ig- 
nitibility. 

Viscous  Heating  Cate:  a  <  0 

We  now  consider  the  case  when  viscous  heating  exceeds  the  ab¬ 
solute  value  of  the  temperanire  difference  between  the  two 
streams,  i.e.,  1^1  <  p.  In  this  case,  the  frozen  temperature  profile 
(IS)  achieves  its  maximum  value  of  «  T.  4-  aV4p  at  die  location 
^  s  (1  -  p/p)/2,  in  the  interior  of  the  mixing  layer.  The  stnicture 
of  the  mixing  layer,  shown  in  Fig.  4c,  consists  of  two  0(1)  frozen 
outer  zones  separated  by  a  narrow  diffusive-reactive  zone  near 
4  B  The  parabolic  form  of  T/near  ^  e  ^  suggests  that  the  appro¬ 
priate  stretched  inner  coordinate  is  Z  «  Vp  ~  ^)/^>  n>ch  that 

/"(I  -  X) »  /"(T1,)/(l  -  X)  +  0(1)  (25) 


which  is  valid  for  0  <  a  <  1 .  Substituting  Eq.  (24)  into  Eq.  (23),  we 
obtain  an  explicit  formula  for 

Equation  (23)  reveals  that  the  leading-order  value  for  the  nondi- 
mensional  ignition  distance  ^  is  a  function  of  p,  which  groups  the 
two  parameters  Af.  and  X.  There  is  no  additional  dependence  on  X 
so  that,  in  the  absence  of  viscous  heating  (Af.  s  0),  is  unaffected 
by  changes  in  X  to  leading  order.  This  implies  that  variations  in  the 
flowfield  at  the  cold  boundary  have  only  a  secondary  effect  (see 
Appendix  B).  in  agreement  with  the  previous  numerical  results. 

We  also  note  that  the  structure  equation  (20)  is  locally  similar. 
Furthermore,  the  matching  condition  (22)  shows  that  there  is  no 
heat  transfer  between  the  inner  structure  and  the  outer  frozen  flow 
up  to  0(e).  Consequently,  the  effect  of  nonsimilarity  is  of  higher 
o^r  and  does  not  appear  in  the  leading-order  analysis,  as  was  true 
for  the  flat-plate  flow.*  However,  the  numerical  result  shown  in 
Fig.  3a  seems  to  suggest  that  there  is  a  larger  difference  between 
the  similar  and  nonsimilar  solutions  for  the  present  study  as  com¬ 
pared  to  the  “subadiabatic  wall  case”  of  the  flat-plate  analysis.' 
This  can  be  explained  by  examining  the  relative  order  of  magni¬ 
tude  of  the  diffusion  and  (nonsimilar)  convection  terms  in  the 
equations  governing  the  structure  for  each  flow  configuration.  It 
can  be  readily  shown  that  the  ratio  of  convective  to  diffusive  terms 
in  the  structure  equation  is  0(l/C.£-2)  for  the  present  problem, 
which  is  larger  than  the  C^e’)  ratio  for  the  flat-plate  problem.  As 
discussed  in  the  previous  section,  the  physical  implication  is  that 
the  streamwise  flow  velocity  at  the  ignition  location  is  relatively 
larger  for  the  present  problem  than  for  the  flat-plate  problem. 


where  [1  - /'(TlflJfd  -  X)  =  ^.  We  seek  a  solution  for  the  temper¬ 
ature  in  the  reaction  zone  of  the  f^rm 

TtofZ.  0  «  r,(Z)  +  eeo(Z.  C)  +  e«ei(Z.  C)  +  (26) 

To  achieve  proper  balancing  of  the  diffusion  and  reaction  terms,' 
we  introduce  the  rescaled  Damkbhler  number  as 


A  = 


1 

M[r(Ti,)/(i-x.)i’ 


exp(-r,/r,) 


(27) 


After  substimting  the  expansions  into  Eq.  (16),  the  structure  equa¬ 
tions  at  the  first  two  ordm  become 


az^ 


-AexpfBo-Zb 


(28) 


As  discussed  in  the  flat-plate  study,'  integration  of  the  preceding 
equations  provides  boundary  conditions  and  a  jump  coridition  for 
the  outer  temperature  profile,  which  is  expanded  in  each  region  as 


C) »  Tf<X)  f  ;)  +  CHe”) 


(29) 
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We  substitute  Eq.  (29)  (16)  and  introduce  the  new  stream- 

wise  coordinate  Ts  r>.(Vnii A)  to  obtain  an  equation  for  the  lead¬ 
ing-order  temperature  perturbation  of  the  form 


3‘Oo  211- (1-X)^]  , 

d4'  ir(T1)/(l->.)]''5^ 


(30) 


Matching  solutions  of  Eq.  (30)  with  the  inner  solutions  (28)  yields 
the  conditions 


0^(0.  T)=«>;(i.T)=«>;(t-“)=o 


(31) 


* -exp(«>o(5,.T) +T]  (32) 

The  system  (30-32).  which  has  also  been  obtained  independently 
by  Grosch  and  Jackson,^  has  a  similar  form  to  that  appearing  in  our 
previous  study.'  The  solution  of  the  present  system  exhibits  ther¬ 
mal  runaway  at  T;,  which  is  determined  as  a  function  of  the  two  pa¬ 
rameters  and  X.  The  numerical  solution  of  t/  as  a  function  of  ^ 
for  several  values  of  X  is  shown  in  Fig.  6.  From  the  definition  of  t/, 
we  can  readily  determine  the  ignition  distance  as 

»  exp[t,(4,,X)] 

We  remark  that  in  the  present  analysis,  nonsimilar  effects  have 
been  accounted  for  in  the  outer,  frozen  regions,  and  thus  the  solu¬ 
tion  possesses  evolutionary-type  behavior  similar  to  that  of  ther¬ 
mal  explosion.  If  the  nonsimilar  term  d/dr  in  Eq.  (30)  is  neglected, 
then  the  equation  can  be  easily  integrated  to  yield  an  explicit  ana¬ 
lytical  solution  which  possesses  a  turning  point  at 


Fig.  7  Functional  relation  between  the  ignition  DarakObler  mmiber 
and  r  for  the  intermediate  cam  for  g  ^  l. 

sources  generate  comparable  amounts  of  heat.  We  now  consider 
this  limiting  case  to  provide  a  complete  description  of  ignition  in 
compressible  mixing  layers. 

We  rescale  a  as  a  s  Jia,  introduce  the  stretched  coordinate  X 
s  Jiil  /Jt,  and  seek  solutions  for  the  inner  temperatuie  profile  as 

rjx.C)  =  r/(X)4-eeo(x.C)  +  0(e“)  (35) 

To  leading  order,  the  inner  structure  equation  is  found  to  be 

-A„X*exp(eo-X’-rX)  (36) 

3X 

where  F  =  a/Jv^.and 

=  ^^yr.X-'’*’*’exp{-r./rj  (37) 

e-£  Vii 


T,(^)=-l-C„[^(l  -^)].  (34) 

This  locally  similar  result,  however,  underestimates  the  ignition 
distance  by  a  substantial  amount  as  shown  by  the  vertical  dotted 
lines  C  and  D  in  Fig.  3b. 

Intermediate  CaM:  |a|«  1 

Neither  of  the  precediM  analyses  is  valid  when  Ial4:  1,  or 
more  specifically  which  occurs  when  the  two  ignition 


5. 


Fig.  <  Functional  relation  between  the  streamw  be  x,  and  transverse 
4r  location  of  ignition  for  the  viscous  beating  case  for  various  L 


The  outer  temperature  profile  is  expanded  as  in  the  hot-stream  case 
such  that  the  boundary  and  matching  conditions  are  given  by 

de„ 

fio(O.C)*0.  ^  (-,(;)  =  0  (38) 

which  are  identical  to  those  in  Eq.  (22).  The  inner  structure  prob¬ 
lem  for  the  present  case  is  locally  similar,  consistent  widi  die  hot- 
stream  case.  Furthermore,  all  variations  of  with  X  and  M.  are 
felt  through  the  single  lumped  parameter  p.  Numerical  solution  of 
Eqs.  (36)  and  (38)  shows  that  for  a  given  value  of  F,  no  solution 
exists  beyond  a  critical  value  which  is  defined  as  the  ignition 
Damkbhler  number.  Solutions  orily  exist  for  q  >  0,  and  we  dwose 
q  w  I  to  be  consistent  with  the  other  cases  considered  eariier.  In 
Fig.  7,  S„j  is  plotted  as  a  function  of  the  parameter  F  and  is  seen 
to  increase  monotonically  with  F.  A  useful  correlation  for  this 
curve  is  given  by 


A-/(F)«fJ  1  + 


1.141  exp(0.547F) 


1  +  0.433  exp(-0.756F- 0.01 7r*). 


(39) 


Given  the  value  of  A„;.  the  ignition  distance  for  the  intermediate 
case  is  determined  from  the  definition  of  A„  in  Eq.  (37)  as 


(40) 


in  which  q  =  1  when  the  correlation  (39)  is  used. 
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SwMMr}'  of  AfynqHotk  Analysis 

Equations  (23),  (33),  and  (40)  now  provide  a  complete  descrip¬ 
tion  of  the  ignition  distance  over  a  wide  range  of  parameters.  We 
letnark  that  we  have  not  treated  the  case  of  c'He)  temperature  and 
velocity  differences  as  studied  by  Jackson  and  Hussaini.^  For  that 
limiting  case,  reaction  occurs  over  the  entire  domain,  and  a  full  nu¬ 
merical  approach  must  be  employed.  We  also  note  that  the  present 
study  can  be  extended  to  describe  a  three-step  chain  branching-ter¬ 
mination  mechanism  proposed  by  Biriian  and  Law.'^  For  such  a 
case,  in  the  limit  of  fast  recombination  in  which  the  radicals  are  in 
steady  state,  the  problem  degenerates  to  that  of  the  one-step  mod- 
el,*-'^and  thus  tli«  present  analysis  can  be  readily  applied. 

The  predicted  ignition  distances  obtained  in  the  analyses  for  the 
inevious  three  cases,  together  with  the  direct  numerical  results,  are 
plotted  in  Fig.  8  as  a  function  of  Af..  To  avoid  cluttering  the  graph, 
only  the  nonsimilar  results  are  included.  The  curves  demonstrate 
how  the  three  different  cases  match  together  to  provide  a  global 
picture  for  the  ignition  distance  as  a  function  of  Mach  number.  The 
overall  agreement  with  the  numerical  results  is  also  seen  to  be 
quite  good.  For  the  case  of  X  s  1,  which  corresponds  to  a  uniform 
parallel  flow,"  there  is  no  viscous  heating,  and  thus  MJii  is  deter¬ 
mined  by  the  result  of  the  hot-stream  case,  Eq.  (23),  for  all  A/..  As 
discuss^  earlier  in  Fig.  S,  the  ignition  distance  in  this  case  varies 
linearly  with  the  freestream  velocity,  as  illustrated  by  the  top  curve 
in  Fig.  8.  For  nonunity  values  of  X,  however,  MJ^  passes  through 
a  maximum  value  and  then  decreases  rapidly  as  a  result  of  viscous 
heating.  Since  the  amount  of  viscous  heating  is  propoitional  to  Af_ 
and  the  reaaion  rate  increases  exponentially  with  the  characteristic 
temperature  of  the  reaction  zone,  we  observe  a  drastic  decrease  of 
the  ignition  distance. 

The  maximum  point  (Af_^)„„  of  each  curve  in  Fig.  8  may  have 
practical  importance  in  the  design  of  the  scramjet  engine,  as  it  rep¬ 
resents  the  most  difficult  situation  in  which  to  achieve  ignition. 
The  value  of  (Af.^)„„  can  be  used  to  evaluate  the  critical  length 
of  a  flame  holding  region  that  is  sufficient  to  ensure  ignition  for  all 
Mach  numbers.  In  holding  regions  of  smaller  dimension,  it  may  be 
necessary  to  vary  operating  conditions  to  achieve  shorter  ignition 
distances.  We  will  now  discuss  how  the  critical  distance  varies 
with  each  of  the  system  parameters. 

Since  this  maximum  point  is  described  by  the  hot-stream  case, 
the  critical  ignition  distance  is  determined  from  Eqs.  (23)  and  (S). 
It  is  clear  from  Eq.  (23)  that  the  actual  distance  is  most  sensitive  to- 
variations  in  the  hot-stream  temperature  T..  An  increase  in  will 
result  in  a  substantial  decrease  in  the  critical  ignition  distance  due 
to  the  Arrhenius  dependence. 

To  assess  the  effect  of  the  two  remaining  relevant  parameters 
and  Af„  we  fix  both  and  X«  1 ).  Then  it  is  convenient  to 
work  with  the  scaled  variable 

yP 

e  ’f«e 


Fig.  •  Globai  log-log  ploi  of  the  predicted  ignition  disunce  as  a  func- 
Mon  of  the  frecstr.’am  .Mach  number  for  various  X. 


P 


Fig.  9  Maximum  ignition  length  in  terms  of  the  nondimcnsioiwl 
quantity  L!  defined  in  Eq.  (41)  as  a  function  of  p. 

which  corresponds  to  the  properly  normalized  distance.  It  follows 
from  Eq.  (23)  that  L,  is  a  function  of  only  two  parameters,  P  and  |i, 
representing  and  Af.,  respectively.  For  q  s  l,  the  correlation 
function  (24)  can  be  used  to  t^tain  an  explicit  expressitm  for  as 

L, (p.  p)  =  2c-2>/p  (p  -  pK2  -  P  +  uVd  -  P  +  p)*  (42) 

When  P  is  fixed,  Lj  attains  a  maximum  at  a  critical  value  of  p, 
which  represents  the  Mach  number  at  the  maximum  poim  of  each 
curve  in  Fig.  8. 

In  Fig.  9  we  have  plotted  the  maximum  value,  L; ,  as  a  function 
of  p.  Since  Q  and  P  both  vary  with  T.,  it  is  necessary  to  hold  f. 
fixed  when  interpreting  Fig.  9.  Therefore,  an  increase  in  p  here  im¬ 
plies  a  decrease  in  the  cold  boundary  temperature  which  leads 
to  an  increase  in  the  heat  loss  from  the  hot  boundary.  Conse¬ 
quently,  ignition  is  more  difficult  to  achieve,  and  L}  is  seen  to  in¬ 
crease  with  p. 

Conclusions 

The  thermal  ignition  in  a  supersonic  mixing  layer  between  two 
streams  of  reactants,  allowing  for  aihitrary  differences  in  both  ve¬ 
locities  and  temperatures,  has  been  studied  numerically  and  as¬ 
ymptotically.  The  asymptotic  analysis  has  revealed  that  die  reac¬ 
tion  zone  structure  takes  different  forms,  depending  on  whether  die 
primary  source  for  ignition  is  the  hot  freestream,  the  viscous  heat¬ 
ing  generated  within  the  shear  layer,  or  a  combination  of  the  two. 

Similar  to  the  flat-plate  case,'  the  ignition  behavior  in  die  mix¬ 
ing  layer  is  characterized  by  the  parameter  a,  which  is  a  measure 
of  the  relative  effects  of  the  temperature  difference  of  die  two 
freestreams  and  the  amount  of  viscous  heating.  In  die  presem  mix¬ 
ing  layer  situation,  however,  there  is  an  additional  parameter  in¬ 
volved.  namely,  X,  which  represents  the  velocity  ratio  of  die  slow¬ 
er  to  the  faster  freestreams.  In  the  subsonic  limit  (A#.  *  0),  viscous 
heating  is  negligible  and  an  increase  in  X  results  in  a  small  increase 
in  the  ignition  distance.  For  finite  values  of  Af^  such  diat  viscous 
heating  is  effective,  a  direct  comparison  of  the  resuitt  for  variotts 
values  of  X  is  difficult  because  the  amount  of  viscous  heating  de¬ 
pends  on  both  Af.  and  X.  As  a  result  of  the  combined  effect  of  the 
two  parameters,  the  ignition  distance  exhibiu  several  different 
types  of  behavior  as  Af.  is  varied.  In  particular,  for  the  hot-stream 
case  (a  >  0).  variations  in  X  produce  only  higher  order  effectt  and 
the  ignition  characteristics  are  structurally  similar  to  those  of  a  uni¬ 
form  parallel  flow."  For  the  case  of  a  <  0,  viscous  heating  is  die 
dominant  energy  source  for  ignition,  and  the  point  of  maximum 
temperature  shifts  to  the  interior  of  the  mixing  layer.  Conse¬ 
quently.  variations  in  the  flowfleld  at  both  boundaries  have  a  sig¬ 
nificant  influence  on  the  ignition  distaiKe. 

The  effects  of  flow  nonsimilarity  have  also  been  usessed  by 
comparing  solutions  of  the  full  system  with  those  found  by  assum- 
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ing  local  similarity  in  the  flowfield.  Nonsimilar  effects  are  most 
important  when  viscous  heating  is  sufficiently  large  so  as  to  cause 
a  temperature  bulge  in  the  interior  of  the  mixing  layer.  In  general, 
the  locally  similar  assumption  has  been  found  to  underestimate  the 
ignition  distance,  in  agreement  with  the  results  of  the  flat-plate 
case.'  It  has  also  been  demonstrated  that  nonsimilar  effects  are 
more  imponani  for  the  mixing  layer  than  for  the  flat-plate  flow,  in¬ 
dicating  the  importance  in  retaining  these  terms  in  studies  involv¬ 
ing  the  ignition  behavior  of  supersonic  mixing  layer  flows. 

Appendix  A:  Ignition  at  the  Siow  Boundary 
When  P  <  0,  ignition  occurs  near  the  slow  boundary.  However, 
it  is  possible  to  redefme  certain  quantities  to  obtain  the  same  lead¬ 
ing  order  structure  problem  as  in  Eqs.  (20-22).  It  is  convenient  to 
first  choose  the  transverse  coordinate  ^  to  be 


where  the  constants  of  integration.  C.  and  are  functions  of  k. 
Combining  the  preceding  formulas  yields 

(B2) 

where 

1-/-  720.  exp(-2*) 

It  follows  from  Eq.  (B3)  thait 

u£,~-z^-Uz-*-U[j2DJ(\-k)\  (B4) 

which,  to  leading  order,  yields 


so  that  ignition  again  occurs  at  ^  =  0.  Following  the  same  proce¬ 
dure  as  in  the  hot-stream  case,  the  identical  structure  equation.  Eq. 
(20),  is  obtained  if  we  redefine  the  variables  as 


(B5) 

Substimting  Eq.  (BS)  into  Eq.  (B2)  provides  the  result 


.r- ' 


A  * 


(A2) 

Tl‘'”’*'exp(-Vr_) 

(A3) 

Thus  if  use  is  made  of  Eq.  (24)  in  Eq.  (A3),  the  ignition  distance  is 
found  to  be 


of J-r. 

5/  =  (“)  y— txp(T/T^)  (A4) 

which  depends  on  an  additional  factor  X. 

For  the  intermediate  case,  lal  c  1.  we  again  use  the  coordinate 
in  Eq.  (Al),  and  by  redefining  the  variables 

e »  a  =  T_-  7.-  p  (A5) 

we  obuin  the  same  structure  equations  (36)  and  (38).  Now  the  re¬ 
duced  Damkbhler  number  is  given  by 

•  p/2-lyrP 

*  - ,  ''-*1  r - (A6) 

Xr»£  Tii 


-222^2- (B6) 

that  is  used  to  derive  Eq.  (20).  We  note  that  the  last  term  in  Eq. 
(B4)  is  a  subdominant  term  and  thus  was  eliminated  in  the  present 
leading  order  analysis.  Consequently,  the  final  structure  equation 
(20)  is  independent  of  the  parameter 
For  the  special  case  of  a  uniform  flowfield,  i.e.,  X  ^  l,  the  solu¬ 
tion  for  the  flowfield  is  simply/*  t)  and  the  factor  [/"/(I  -  X)]* 
multiplying  the  diffusion  term  in  Eq.  (16)  has  the  explicit  form 

^exp(-2q^)  (B7) 

It  is  convenient  to  woiX  in  terms  of  the  original  coordinate  ^  *  eifc 
(ti/72  )/2  (see  Eq.  (13)]  so  that  the  asymptotic  behavior  of  Eq. 
(B7)  near  the  ignition  point  qj  is  given  by 

^exp(-2Ti') (B8) 

wliere  the  large  number  qj  can  be  evaluated  by  taking  the  loga¬ 
rithm  of  Eq.  (B8),  yielding 

q*-2t«(a/2cx7n)  -  2t«qi  -  - 


so  that,  when  use  is  made  of  the  correlation  (39)  for  A„j,  the  ex- 
pressitm  for  the  ignition  distance  becomes 


(B9) 


^  I.-' — 


(A7) 


in  terms  of  the  inner  variable  x  defined  in  die  hoi-stitam  case. 
When  only  the  leading-order  term  on  the  right-hand  side  of  Eq. 
(B9)  is  retained,  namely,  the  first  term,  we  obtain 


Appendix  B:  Asymptotic  Behavior  of /"/(I  -  X)  as  ^  0 

The  aiuilysis  of  the  hot-stream  case  requires  the  asymptotic  be¬ 
havior  of  /”/(!  -  X)  as  ^  0  which,  from  definition  (13).  corre¬ 

sponds  to  q  ->  M.  By  integrating  Eqs.  (1 )  and  (7)  near  q  *  <».  it  oin 
be  shown  that’ 

/’'-O.exp(-r’)  (Bl) 


qj--f.62  (BIO) 

which  is  consistent  with  our  previous  result  Eq.  (B6),  whereas 
Lihfn  and  Crespo' '  retained  the  second  term  on  the  right-hand  side 
of  Eq.  (B9).  A 'though  the  two  expressions  must  be  asymptotically 
the  same  in  the  limit  e  ->  0,  in  practice  they  can  result  in  quite  a 
large  difference  in  the  predict^  ignition  distance  as  sht^  in 
Fig.  3a. 
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